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We investigate the nature of the Mott-insulating phases of half-filled 2Ai'-component fermionic cold atoms 
loaded into a one-dimensional optical lattice. By means of conformal field theory techniques and large-scale 
DMRG calculations, we show that the phase diagram strongly depends on the parity of A'^. First, we single out 
charged, spin-singlet, degrees of freedom, that carry a pseudo-spin S — N/2 allowing to formulate a Haldane 
conjecture: for attractive interactions, we establish the emergence of Haldane insulating phases when A*' is even, 
whereas a metallic behavior is found when A'^ is odd. We point out that the TV = 1, 2 cases do not have the 
generic properties of each family. The metallic phase for A*' odd and larger than 1 has a quasi-long range singlet 
pairing ordering with an interesting edge-state structure. Moreover, the properties of the Haldane insulating 
phases with even TV further depend on the parity of TV/ 2. In this respect, within the low-energy approach, we 
argue that the Haldane phases with TV/ 2 even are not topologically protected but equivalent to a topologically 
trivial insulating phase and thus confirm the recent conjecture put forward by Pollmann et al. [Pollmann et al, 
[MXiv:0909.4059| (2009)]. 
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I. INTRODUCTION 

Topological phases have attracted much interest in recent 
years due to their robustness against perturbations and their 
relevance to quantum computation. A topological ordered 
phase is a gapped phase which displays a protected ground- 
state degeneracy dependent on the topology of the manifold 
in which the model is embedded-l^ This phase is not charac- 
terized by a local order parameter and falls beyond the usual 
symmetry breaking paradigm of condensed matter physics.^ 

One of the simplest examples of topologically ordered 
phases is the Haldane phase in quantum spin chains. In 1983, 
Haldane argued that the spin-S* Heisenberg chain displays 
striking different properties depending on the parity of 25'.'^ 
While half-integer Heisenberg spin chains have a gapless be- 
havior, a finite gap from the singlet ground-state (GS) to the 
first triplet excited states is found when 25 is even. On top 
of the existence a gap, the spin-1 phase (the so-called Hal- 
dane phase) has remarkable exotic properties which may be 
regarded as manifestations of the existence of a topological or- 
dered phase. This phase is not characterized by a local order 
but displays non-local string long-range ordering which sig- 
nals the presence of a hidden Neel antiferromagnetic order."* 
The latter can be revealed through a non-local unitary trans- 
formation and the emergence of a complete breaking of a Z2 
X Z2 symmetry.^ One remarkable resulting consequence of 
the Haldane phase is the liberation of fractional spin- 1/2 de- 
grees of freedom at the edge of the sample when the chain is 
doped by non-magnetic impurities.l^ 

Haldane's conjecture is now well understood and has been 
confirmed experimentally in quasi- ID compounds as well as 
numerically (see for instance Refs,7 8). The Haldane phase 



displays unusual and interesting physical properties so that 
it is important to experimentally stabilize it in other con- 
texts. In this respect, it has been argued that the Haldane 
phase is relevant to Josephson junction array systems.^ Fur- 
thermore, it is likely that the Haldane physics will be ex- 
plored experimentally in the near future in trapped ultracold 
atomic systems thanks to the tunability of interactions in these 
systems using optical lattices and Feshbach resonances. A 
first possible direction is to consider spin-1 bosons loaded 
into a one-dimensional (ID) optical lattice with one atom 
per site so that the Haldane phase is one of the possible in- 
sulating phases of this model. "^ A second route consists in 
preparing ID ultracold quantum gases with dipolar interac- 
tions, like ^^Cr bosonic atoms, where a Haldane insulating 
(HI) phase has been predicted.'' '"* Finally, we have recently 
shown that a similar phase can also be stabilized by consid- 
ering ID spin-3/2 cold fermions at half-filling with contact 
interactions onlyP 

In this paper, we pursue our investigation of the HI phase 
in the context of ID ultracold fermionic alkaline atoms in 
the general half-integer (hyperfine) spin F ~ N — 1/2 case 
at half-filling (TV atoms per site). In this respect, we will 
use complementary analytical (renormalization group (RG) 
analysis, conformal field theory (CFT))'* and density-matrix 
renormalization group (DMRG)'^ techniques to fully deter- 
mine the nature of the Mott-insulating phases at half-filling 
when TV > 2. The starting point of the analysis is the lattice 
model of 2TV components cold fermions with contact inter- 
actions. Due to Pauli principle, low-energy s-wave scattering 
processes of spin-F fermionic atoms are allowed only in the 
even total spin J = 0, 2, ... , 2TV — 2 channels, so that the 
general effective Hamiltonian with contact interactions reads 



as follows in absence of a magnetic fieldP^ 
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where c]^ ^ is the fermion creation operator corresponding to 
the 2N hyperfine states (a = 1, . . . , 2N) at the i* site of the 
optical lattice. The pairing operators in Eq. (fill are defined 
through the Clebsch-Gordan coefficients for spin-i^ fermions: 
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In the general spin- 



F case, there are N coupling constants Uj in model (fTJ, 
which are related to the N two-body scattering lengths of 
the problem. In the following, in order to simplify the anal- 
ysis of the Mott-insulating phases when iV > 2, we per- 
form a fine-tuning of the different scattering lengths in chan- 
nels J > 2, i.e. U2 = ••• = t/2Ar_2_. Using the identity 
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with U = 2U2, V ^ Uo - U2, and m 
Sq cL iC^ i is the density at site i. In Eq. 
BCS pairing operator for spin-F fermions is 



J2|l, the singlet 



Ea/3 cl,. Ja/j4^^ = - Ea (-1)" cl^,4w+i-a,»' the matrix 
J being a 2N x 2N antisymmetric matrix with J^ = — /. 



When iV = 1, Pqq ^ coincides with the Cooper pairing c. 






SO that model (|2]i is equivalent to the spin- 1/2 Hubbard model. 
Model (|2]l obviously conserves the total number of fermions 
- no atoms are dynamically created or annihilated. This 
conservation law is associated to a U(l) continuous sym- 
metry Ca,i ^" e^^Ca,i, which, by analogy with condensed 
matter dealing with electrons, we will refer to as a "U(l)c 
charge symmetry", the charge being simply the number of 
fermions. On top of this symmetry, model (|2]i displays an 
extended continuous symmetry for A^ > 1 in spin space. 
When V = 0(Uo = U2) model Q is the Hubbai'd model for 
2iV-component fermions with a U(2iV)= U(l)c x SU(2A^) 
invariance. The Hamiltonian (|2]i for V ^ Q still displays 
an extended symmetry since the BCS singlet-pairing operator 
Pqq j is invariant under the Sp(2A^) group which consists of 
27V X 2N unitary matrices U that satisfy U*jW = J. When 
V^ 7^ 0, the co ntinuous symmetry of model (J2]i is thus U(l)c x 
Sp(2N)^^^ln the F = 3/2 case, i.e. iV = 2, there is no fine- 
tuning; models (fill and (J2]i are equivalent and share an exact 
Sp(4)~ SO(5) spin symmetry.^' The zero-temperature phase 
diagram of model (J2]) away from half-filling has been inves- 
tigated by means of a low-energy approach^^ ^'^ in the gen- 
eral N case, and by Quantum Monte Carlo and DMRG cal- 
culations for N — 2}^ ^^ A rich exotic physics emerge when 



N > 2 with, in particular, the stabilization of a superconduct- 
ing instability with charge 2Ne for attractive interactions and 
at sufficiently low density. ^^^"^ 

At half-filUng (when fi = NU + V/N), model ^ enjoys 
a particle-hole symmetry c^^i — > i^^Y ^r J^apct ^ which 
plays a crucial role in the following. In the N — 1 case, 
it is well known that the particle-hole symmetry enlarges the 
U(l)c charge symmetry of the s piri-1 /2 Hubbard model to an 
SU(2)c symmetry at half-fiUing.E^ED in addition, the physics 
of the half-filled spin- 1/2 Hubbard model for repulsive and 
attractive interactions are related through a canonical trans- 
formation c^i -^ (— l)*c| j, c^^i — >■ c^.j. While for U > 
a Mott-insulating phase with one gapless spin modes is sta- 
bilized, there is a spin gap for attractive intera ction which 
marks the emergence of a singlet-pairing phase. l^l^ When 
A^ > 1, all these properties do not generalize; in particular, 
the symmetry enlargement of the charge degrees of freedom 
at half-filling requires an additional fine tuning V = NU to 
display an SU(2)c x Sp(2A^) global invariancclIS We have 
shown in Ref IT5l that this SU(2)c symmetry is central to the 
emergence of an even-odd scenario for attractive interactions 
in close parallel to the famous Haldane conjecture in spin- 
S SU(2) Heisenberg chains. In this respect, we have identi- 
fied a spin-singlet pseudo-spin N/2 operator which governs 
the low-energy properties of the model in the vicinity of the 
SU(2)c line for attractive interactions. This operator gives rise 
to a Haldane-charge conjecture with the emergence of a HI 
phase when N is even, while a metallic phase is stabilized 
when N is odd. Such a scenario has been checked in ReffTSJ 
by a low-energy approach in the N — 2 case and DMRG cal- 
culations for A^ = 2, 3 in the vicinity of the SU(2)c line. In the 
special N = 2 case, these complementary techniques reveals 
unambiguously the existence of a HI phase with non-local 
string charge correlations and pseudo-spin- 1/2 edge states. 

In this paper, we extend the results of our letter Ref 15 by 
determining the zero temperature phase diagram of model (J2]) 
at half-filling by means of a low-energy approach in the gen- 
eral N case and DMRG calculations for TV = 2,3,4. On 
top of the confirmation of the Haldane-charge conjecture, we 
show that the iV = 1 and N — 2 cases are special and are not 
the generic cases of each family. In particular, f or A^ > 1 odd, 
the metallic phase with dominant singlet-pairing correlation 
has an interesting edge-state structure when open-boundary 
conditions (OB C) are used, similarly to the spin-3/2 Heisen- 
berg chain.EiEIIpQj- ^n A^ even > 2, a new gapless phase with 
dominant singlet-pairing instability is stabilized between the 
HI phase and the rung-singlet (RS) phase. In addition, we 
show, within the low-energy approach, that the HI phase has 
striking different properties depending on the parity of A^/2. 
When A^/2 is even, the HI phase turns out to be equivalent 
to the topologically trivial RS insulating phase whereas it is 
a topologically ordered phase when N/2 is odd in full agree- 
ment with the recent findings in the study of integer Heisen- 
berg spin chains.-'^ ^* 

The rest of the paper is organized as follows. In Sec. II, 
we discuss the strong-coupling analysis of model (|2]) along 
special highly-symmetric lines which give some clues about 
the nature of the Mott-insulating phases. The low-energy ap- 



proach of the general N case is presented in Sec. III. In Sec. 
IV, we map out the phase diagram of model (|2]i with N — 2 
by means of intensive DMRG calculations. Section V and 
VI describe our DMRG results respectively for the iV = 3, 4 
cases to complement the low-energy approach. Finally, our 
concluding remarks are given in Sec. VII. 



The energy level {n, k) is ^{n, fc)-fold degenerate, with: 

_ 2iN-n + 2k + l)i2N+l)l 
^''~in^2ky.{2N-n + 2k + 2)\- ^' 

At half-filling, n is set by the particle-hole symmetry: /i = 
NU + V/N, and the energy levels read: 



II. STRONG-COUPLING ANALYSIS 

Before investigating the zero-temperature phase diagram 
of model (J2|i by means of the low-energy and DMRG 
approaches, a strong-coupling analysis along the highly- 
symmetric lines of the model might be useful to shed light 
on the possible Mott-insulating phases. To this end, let us first 
consider the energy-spectrum for the single-site problem. 

The Hubbard term of Eq. (J2]) distributes the different states 
into energy levels with the same number of particles n, with 
n = 0, . . . , 2N; this is the one-site spectrum of the U(2A^) 
Hubbard model. The singlet-pairing term in Eq. (J2]) with cou- 
pling constant V will split these levels into levels with differ- 
ent pairing schemes, denoted {n,k). The level {n,k) group 
states with n particles, among which 2fc particles are in k 
Sp(2A^) singlets. These states transform in the ujn-2k rep- 
resentation of Sp(2A^). Note that, for a given number of parti- 
cles n, < fc < E{n/2) if n < N,iindn-N < k < E{n/2) 
if n > N, where E{x) is the floor function. In order to write 
down the eigenstates in terms of fermionic operators, let us 
define the pair operator that creates a pair of fermions with 
spins a and 2N + 1 — a, by: 
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of these operators, the singlet pairing operator PJq ^ 
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We now need to define a set of iV — 1 linear combinations 



that we label as IIJ ■ (with I — 



of pIj, "orthogonal" to fJq , 

1, • • • ', iV - 1). The n - 2k particles that are not Sp(2A^) 
singlets then divide into two kinds: they can be either written 
as linear combinations of pairs of particles with spin {a, 2N + 
1 — a) and thus as a combination of pair operators HJ ^, or they 
are unpaired and can be only written with a single creation 
operator c]^ ^. In the end, the eigenstates that belong to the 
energy level (n, fc) are written as: 
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(5) 
where m labels the state, Ain. k.m is a normalization factor, 
p is the number of "single" particles, and 2q is the number of 
"paired" particle that cannot be penned down in terms of pJq ^. 
The energy of the eigenstates (jsl) only depends on (n, k) and 
reads: 
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At this point, we can consider two important highly- 
symmetric lines for all N: V = (respectively V = 
NU) with the emergence of a U(2iV) (respectively SU(2)c x 
Sp(2A^)) extended symmetry. We also mention that in the spe- 
cial N — 2 case, there is an additional SO(7) symmetric line at 
half-filling when V = — 2 t/.E2 However, despite the fact that 
we indeed find an additional degeneracy for the one-site prob- 
lem in the general N case on the special line V — —N'^U/2, 
the latter seems not to correspond to an enlarged symmetry 
since the kinetic term lifts it. 



A. Strong-coupling argument close to the V = line. 

When F = 0, as already stated in the introduction, model 
Q is equivalent to the U(2iV) Hubbard model. The degen- 
eracies of the energy-spectrum dSll with V^ = are related to 
the dimensions of representations of the SU(2A^) group. In 
particular, when [/ > 0, we observe from Eq. ([8]) that the 
lowest-energy states correspond to n = A^ and transform in 
the antisymmetric self-conjugate representations of SU(2A^) 
(representation described by a Young tableau with one col- 
umn of N boxes). This case has been studied in Ref. [35] 
and in the strong coupling limit the model is equivalent to 
an SU(2A^) Heisenberg spin chain where the spin operators 
belong to the antisymmetric self-conjugate representation of 
SU(2A^). The latter model is expected to have a dimerized 
or Spin-Peierls (SP) two-fold degenerate GS, where dimers 
are formed between two neighboring sites.EESlin the N = 2 
(i.e. SU(4)), such a SP phase has been ascertained by means 
of a low-energy approach. Quantum Monte Carlo and DMRG 
calculations.-'^ "''^ The strong U limit gives the opportunity to 
get a simple physical picture of the GS as well as of the low- 
lying excitations. The two-fold degenerate GS allows for kink 
configurations that interpolate between the two vacua. As 
depicted in FigfTl these kinks have zero charge but carry a 
non-zero SU(2A^) spin since they transform in the antisym- 
metric self-conjugate representation of SU(2A^). Note that 
the system also allows for charged kinks, that carry charge 
Qk = k — N with k = 0, . . . , 2A^. These states trans- 
form in Wfc, the antisymmetric representation of SU(2iV) with 
Young tableau made of a single column with k boxes. Al- 
though at large U they are expected to have a large gap of 
order A^j ~ U {N — fc)^/2 as seen from (Isjl, we nevertheless 
introduce them here since they will play an important role at 



small U (see Sec. [IIIQ . Notice also that these kink excita- 
tions have a collective nature, i.e. their quantum number can- 
not be reproduced by states built by using a finite number of 
fermions. 



iV fermions 
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FIG. 1: Sketch of the kinks supported by the U(2A'^) repulsive Hub- 
bard model, that interpolate between the two degenerate dimerized 
vacua (the boxes indicate an SU(2A'') singlet made of 2N fermions). 
Here A'^ = 5. At large U , in the low-energy sector, all sites have 
exactly TV fermions. (a) Neutral kinks, that are the only low-energy 
excitations at large [/, and transform in the antisymmetric representa- 
tion ujn. Note that they are their own antiparticle. (b) Charged kinks, 
that play a role at small U, depicted here with charge Q = — 2. 

In the attractive case ([/ < 0), the lowest energy states are 
the empty and the fully occupied state, which is an SU(2A^) 
(and Sp(2A^) as well) singlet. At second order of perturbation 
theory, the effective model is thus*^ 



Heff = 
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The first term introduces an effective repulsion interaction be- 
tween nearest neighbor sites. This leads to a fully-gapped 
charge-density wave (CDW) where empty and fully occupied 
states alternate. This phase has a long-range order and is two- 
fold degenerate. 



B. Strong-coupling argument close to the V = NU line. 

The second highly-symmetric line corresponds to the addi- 
tional SU(2)c symmetry in the charge sector for V — NU that 
we have identified in Ref 15) On this line, one easily verifies 
from Eq. i8\, that all pure (pJq j)'^ states (i.e., the states with 
n = 2k, k = 0, . . . , N) are degenerate, with energy E ^ 0. 
Let us give the proper normalization factor for these states: 
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They transform in the spm-N/2 representation of SU(2) and 
we define the coiTesponding pseudo-spin operator acting on 
them as: 



S! = in,-N)/2. 



(11) 



This operator carries charge and is a Sp(2iV) spin-singlet. It 
generalizes the ry-pairing operator introduced by Yang for the 



half-filled spin-1/2 (i.e. N = 1) Hubbard modePorby An- 
derson in his study of the BCS superconductivity. E2 It is easy 
to observe that Si satisfies the SU(2) commutation relations; 
Sj allows to construct the whole set of states with E = from 
|0) with: 



SJ\Po%,) = ViN-k){k + i)\p!;+l) 
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Let us check the commutation relation of S with the Hamil- 
tonian. For the interacting part alone ("Hint = H [< = 0]), we 
have (for a generic filling): 
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(13) 



so that they commute only at half-filling and when V = NU\ 
as for the hopping term, Ut = -^Ej,q (cLc„,i+i + H.c.j, 
it commutes with the total charge pseudo-spin operator if we 
define it as: 
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(14) 



The pseudo-spin operator thus generates a higher SU(2)c 
X Sp(2Af) symmetry at half-filling along the line V — NU 
and we can recast the interacting (on-site) Hamiltonian as 
T^int = 2L/^J5^2 _ ^(^ j^ 2)/4]; the pseudo-spin Si is a 
spin-iV/2 operator. The existence of such an extended SU(2) 
symmetry in the charge sector for N ^ 2 has been first no- 
ticed in Ref. 42 

For a strong attractive U, one can derive an effective Hamil- 
tonian in the strong-coupling regime |f/| ^ t using the stan- 
dard strategy. ^^^ To second order of perturbation theory, one 
obtains the effective model: 
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On the SU(2)c symmetric Hne (V = NU), D = 0, and 



model (15 1 is the spin-A^/2 antifeiTomagnetic SU(2) Heisen- 



berg chain. From this strong-coupling approach, we thus ex- 
pect the emergence of an even-odd dichotomy for attractive 
interactions along the SU(2) line. For even N, i.e. integer 
pseudo spin, the HI phase is formed while a metallic phase is 



stabilized when N is odd, i.e. half-integer pseudo spin. This 
is the same as Haldane's conjecture for model ( [T5] l except that 
the underlying spin S is non-magnetic and carries charge. In 
this respect, we coin it Haldane-charge conjecture. When we 
deviate from this SU(2)c line, the SU(2)c charge symmetry is 
broken down to U(l)c and the single-ion anisotropy appears. 
The phase diagram of the resulting model for general N is 
known from the bosonization work of Schulz.'*'* For even N, 
on top of the Haldane phase, Neel and large-D singlet gap- 
ful phases appear. Using the expression of the pseudo-spin 



operator (111, the Neel and large-D singlet phases correspond 
respectively to the CDW and RS phases. When N is odd, gap- 
less (XY) and gapful (Ising) phases are stabilized in the vicin- 
ity of the SU(2) line. The gapless XY phase can be viewed as 
singlet-pairing phase since SJ ^ Pjg ^. 



III. LOW-ENERGY APPROACH 

In this section, we present the low-energy description of the 
model (J2]) in the general A^ case. This will lead us to map out 
the phase diagram at zero temperature and to show the emer- 
gence of the Haldane-charge conjecture in the weak-coupling 
limit. As it will be shown, the N — 2 case, which was already 
presented in Ref 15, turns out to be very particular and is not 
the generic case of the even N family. 



A. Continuum limit 

The low-energy effective field theory of the lattice model 
(J2]i is derived by taking the standard continuum limit of the 
lattice fermionic operators c^ i, writt en in terms of left- and 
right-moving La , Ra Dirac fermionsP^Eni 
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with X = iao (aq being the lattice spacing) and kp = 7r/(2ao) 
is the Fermi momentum. In the continuum limit, the non- 
interacting part of the Hamiltonian Q corresponds to the 
Hamiltonian density of 2 A^ free relativistic massless fermions: 



Ho = -ivF {Rida;Ra " L^d^L^ 



(18) 



where vp — 2too is the Fermi velocity and we assume in the 
following a summation over repeated indices. The continuous 
symmetry of the non-interacting part of the model is enlarged 
to SO(47V)|l xSO(47V)|fl, since 2iV complex (Dirac) fermions 
are equivalent to 4iV real (Majorana) fermions. This SO(4A^) 
symmetry is the maximal continuous symmetry of 2iV Dirac 
fermions. The corresponding CFT is the SO(4A^)i with cen- 
tral charge c = 2A^."' 

The crucial point is now to find a good basis describing 
the low-energy properties of the model. Some simple con- 
siderations on its symmetries guide us to choose the relevant 
conformal embedding of the problem. No spin-charge sep- 
aration at half-filling is expected for iV > 1, and since the 
global symmetry invariance of model (J2]) is Sp(2A^), we need 



to understand how the non-interacting conformal symmetry 
SO(4iV)i decomposes into Sp(2iV)i CFT. The general list of 
conformal embeddings can be found in R ef J45] and the one 
which is directly relevant to our problem is™ 



SO(47V)i - SU(2)jv X Sp(2]V)i 



(19) 



where the SU(2)jv (respectively Sp(2iV)i) CFT has central 
charge c = 3iV/(iV + 2) (respectively c = A^ (27V + 1)/(7V + 

2)). 

The next step of the approach is to express the 2N{4N — 1) 
SO(4iV)i currents, which are made from all Dirac fermionic 
bilinears, in terms of the currents of the SU(2)jv and Sp(2iV)i 
CFTs, in order to write down the effective interacting Hamil- 
tonian in the new basis. To this end, let us consider the follow- 
ing left currents which appear in the lo w-energy description of 
the model away from half-fillingP^El 






LlT^ljLp, the SU(2iV)i spin currents, 
L'lT^pLjs, the Sp(27V)i spin currents, 
LiKpLp, the SU(2iV)i/Sp(27V)i currents, (20) 



where T^ are the iN^ ~ 1 SU(2iV) generators in the fun- 
damental representation, T° being the N(2N + 1) Sp(2A^) 
generators in the fundamental representation, and T' are the 
2iV^ - A^ - 1 remaining generators of SU(2Af). All these 
generators are normalized in such a way that they satisfy 
Tr(r^r^) = S^^/2. We need then to introduce the SU(2)7v 
currents to complete the conformal embedding. In this re- 
spect, one may use the recognition of the SU(2) pseudo-spin 



operator ( 11 1 to consider the following left current: 
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(21) 



where : : stands for the normal ordering with respect to the 
Fermi sea of the non-interacting theory. Note the unusual def- 
inition of the charge current J'l; the 1/2 factor in Eq. (21 1 is 
there to realize the SU(2)Ar algebra?^ 



J£iz)Jt{0) 



N6' 



ah 



^abc 



2-KZ 



Jl (0) 



(22) 



with a,b ~ 1,2,3 and z — vpT + ix (r being the imaginary 
time). 

At this point, we have only defined AN'^ + 2 currents, and 
we now have to introduce the 4A^^ — 2A^ — 2 other pseudo- 
currents in order to take into account the umklapp terms of 
the form K^^kI, KaKp (K — L, R). In this respect, let us 
consider 






t 

a. ^ 0.(3-'-^ (3') 



(23) 



where the generators T^o are such that, together with J^^, 
they form the set of antisymmetric generators of SU(2A^); 
there are N{2N - 1) - 1 of them, so that all the left 2N{AN - 
1) SO(4A^)i currents are described by: Jl,Jl,J^,Jl, J]^, 
with a = 1, . . . , N{2N + 1), and i = 1, . . . , 2N'^ - N - 1. 



With these currents at hands, we can now derive the low- 
energy effective expression of the interacting part of model (|2]i 
at half-filling. The interacting part of this low-energy Hamil- 
tonian can then be deduced by symmetry, simply by requiring 
the Sp(2iV) invariance: 



Hi, 



54 
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(^^ 



92jyi 

"+H.c.) + 



gzJkJl 



H.c.),(24) 



where we have neglected four-fermion chiral interactions 
which would only introduce a velocity anisotropy. A direct 
continuum limit leads to the identification: 

53 = |ao (u(2N - 1) 

54 = 2[/ao 

55 = ^ao(C/ + 2y). (25) 

Along the U(2A^) line with V^ = 0, we have g\ ==02 = 
^ " - - - iVgs.Wheny = A^[/,modelJ24ldis- 



2y 



-54 



X Sp(2iV) invariance with 92 = 54 



2Ar-l53 

plays a manifest SU(2), 
and 93 — g^. The interaction of model (24i is marginal, so 
that a one-loop RG analysis can be performed to deduce its 
low-energy properties. 



B. Phase diagram for A'^ = 2 

The RG analysis for A^ = 2 has been presented in details in 
Refs. I15I39I For completeness and especially for the compar- 
ison with the DMRG calculations of Sec. IV, we give here a 
brief account of the main results of this case. 

Since Sp(4)i - SO(5)i and SU(2)2 - SO(3)i, the currents 
in Eq. (J24| admit a free-field representation in terms of eight 
Majorana fermions: Cfl li "^ = 1, . . . , 5 describe the fluctu- 
ations of the Sp(4) spin degrees of freedom whereas ^^' £ 
account for the remaining ones. The interacting Hamiltonian 



(24 1 reads as follows in terms of these Majorana fermions: 



'Hint — 
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■52^^?i> ^sc2 
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^^i) 



(26) 



(c^el+^^ii) 54E?^^2+55C^el 



a=l 



One particularity of the Majorana fermion basis is that it al- 
lows for a very simple representation of non-perturbative hid- 
den duality symmetries in the low-energy effective Hamilto- 



nian (26 1. These discrete symmetries are very useful to deter- 
mine the zero-temperature phase diagram.'*^ For model ( 26 1, 




FIG. 2: Phase diagram obtained by the low-energy approach in the 
TV = 2 case. 



we can define two independent duality symmetries: 



^2 



e,' 






?L ~^ ~?Lj 



(27) 



while the right-moving fermions remain invariant. The trans- 



formations (pTji are exact symmetries of Eq. (26 1 if the cou- 
plings are simultaneously changed according to 54.5 — > —54.5 
for r^i, and 52,5 —> —92.5 for the second duality ^2- These 
duality symmetries along with the one-loop RG equations en- 
able us to map out the phase diagram of the N = 2 case. 
This analysis has been done in Refs. I15l39l and we find four 
insulating phases in the phase diagram, depicted in Fig. [2]. 

A first two-fold degenerate phase, which contains the SU(4) 
line with repulsive U, is a SP phase with a non-zero order 

parameter Osp = J2ai~^y(i,,i+i^a-i- Th^ duality sym- 
metry J7i gives a second gapful two-fold degenerate phase 
which is a CDW phase with order parameter Ocdw — 
^^(— l)'cjj, jCa i. This phase contains the SU(4) line with 
negative U, in full agreement with the numerical result of 
Ref. |40l On top of these two-fold degenerate phases, there 
are two non-degenerate insulating phases which are stabilized 
with help of the il2 duality symmetry. Starting from the CDW 
phase and applying the fl2 transformation, one obtains a HI 
phase which includes the SU(2)c line V = 2U with attrac- 
tive U. This gapful non-degenerate phase is equivalent to the 
Haldane phase of the spin-1 Heisenberg chain and displays 
a hidden ordering which can be revealed through a non-local 
string order parameter. This order parameter is built from the 



pseudo-spin operator (111 and the HI phase is characterized 



by the long-range ordering: 



lim 

K-jI-x 



(5fe"^'^-+i'^^-5|) 7^0. 



(28) 



As a consequence of this ordering, this phase displays pseudo- 
spin-1/2 edge states which carry charge but are spin-singlet 
states (holon edge states). '"" Finally, the last insulating phase 
is obtained from this HI phase by applying the ili duality. 
One obtains a gapful non-degenerate RS phase, equivalent to 
the RS phase of the two-leg spin ladder with antiferromagnetic 
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FIG. 3: Ratio of the coupling constants, to wliich tlie RG flow leads 
in the IR limit, for N — 3;6 spans a ring with fixed radius R — O.lt 
in the phase diagram: U — R cos 9, V — R sin 9. Three different 
regions can be defined as function of 9. 



interchain coupling.'22l22l This RS phase has no edge states and 
is characterized by the string-orderingP^l 



lim 



^^-Ei: 



^^•)^0. 



(29) 



Finally, the different quantum phase transitions of Fig. [2] 
can also be determined by means of the duality symmetries 
pTJ l. The transitions SP/CDW and HI/RS are Berezenskii- 
Kosterlitz-Thouless (BKT) transitions with central charge c = 
1 whereas SP/RS and CDW/HI transitions belong to the 2D 
Ising universality class with central charge c = 1/2. 



C. Renormallzatlon Group analysis - General N case 

We now turn to the general N > 2 case which is much more 
involved. The leading effects of the current-current interaction 
of model (|24]i can be inferred from a one-loop RG approach. 
In this respect, it is useful to rescale the coupling constants as: 
51,2.4 —> 27rwF-^^5i,2,4 and 53^5 — ?► 27ruFg3.5, to obtain the 



one-loop RG equations: 

N 



5i 



[iN+l)gl + {N-l)gl + 2{N-l)gl\ 



52 = ^V52 + (^' - iV - 2)g2 + 25455 

53 = i2N'-N-l)gl+gl 

54 = N'^9i9a + {N'^ - N - 2)^254 + .9354 + 5255 



55 



{2N^-N-l)g2g4+g395, 



(30) 



where c/i = dgi/dl{i — 1, ... ,5), I being the RG time. 

Here, we remark the particular character of the case N — 2, 
where the term in N'^ — N ~ 2 cancel out in the equation for 
92 and 54. It is this very term that makes the general N case 
tricky. Indeed, because of it, the duality symmetry il2 of the 



N — 2 case ( 27 1 disappears f or iV > 2 and, with it, fades out 



a very satisfactory way to precisely identify and characterize 
the different phases in the phase diagram. The only duality 
symmetry 51 1 which remains when N > 2, corresponds to the 
transformation La — > iLa on the left-moving Dirac fermions, 
so that 



^1 • ■^r, ^ ^^h ' J). 






(31) 



while the other currents are invariant. Thus, the duality 
transformation 51 1 is an exact symmetry of model ^A\ when 
.94,5 ~^ ^54,5- In particular, the RG Eqs. (30i are indeed 
symmetric with respect to 54.5 — )■ —54,5. 

We have solved numerically the RG equations by a standard 
Runge-Kutta method. In order to have a picture of the RG 
flow, it is useful to draw diagrams on a ring defined by E? — 
C/2 + f2, [/ = RcosO, V = i?,sin6' with6l = . . . 27r; in 
our numerical calculations, we set R = O.lt. The procedure 
is the following: we initiate the couplings gi for a given value 



of 9 and run the Runge-Kutta algorithm on the RG Eqs. ( 30 1. 
The coupling constants gi flow to the strong coupling regime 
under RG time so we need to stop the procedure at one point. 
To this end, we stop the RG iterations as soon as one of the 
coupling reaches a limit value G; this happens at RG time ^0, 
and defines a mass scale A = Oq ^ e^'" that gives an estimate 
of the largest gap in the model (cq ^ is a UV cutoff)- At this 
point, we extract the values of all the gi{lo) and draw them, 
renormalized by 51(^0) (5i/5i(^o)) on the ring diagram. We 
re-initiate the couplings for a new 6 and restart the procedure 
for all values of 6. The resulting diagram looks very similar 
for all A^ > 2 and different regimes can be defined (see Fig. [3] 
for A^ = 3) as a function of the lattice coupling constants U 
and V. Two qualitatively different behaviors of the flow can 
be identified in the asymptotic limit of weak coupling, where 
^0 is large (small Ra^). In regions (I) and (II) of Fig. [3] at large 
Iq, the ratios 5^/51 do not evolve anymore with the RG time 
and have already reached fixed values when the RG iterations 
are stopped. On the other hand, in region (III), gi is always 
the first coupling to reach the limit value G at which we stop 
the RG ; the ratios 5^/51 (^0) for i = 3, 4, 5 vanish in the weak 
coupling limit: 



]im^[gi/giila)] 



0, 



(32) 



whereas the ratio 52/51(^0) remains finite. This property will 
be important when we will derive effective models to describe 



8 



this last region. We now turn to the description of the physical 
properties of the different regimes. 

In the region (I) of Fig. [3] all coupling constants of the low- 
energy effective Hamiltonian (24i flow to strong coupling in 
the infrared (IR) limit at fixed ratio: gi/gi = 1 (z = 1, . . . , 5). 
Along this special isotropic ray, model ([24| displays an ex- 



tended global SO(47V) symmetry and becomes equivalent to 
the SO(47V) Gross-Neveu (GN) model/^ in the sense that the 



low-energy spectrum model of ( 24 1 is adiabatically connected 
to that of the SO(4iV) GN model - for a precise discussion 
see Ref. |48] This phenomenon is an example of a dynamical 
symmetry enlargement by the interactions as found in half- 
filled two-leg Hubbard ladders or in the half-filled U(4) Hub- 
bard chain, with the emergence of an SO(8) symmetry t hat be- 
comes asymptotically exact in the weak-coupling limit.l^^lsoE] 
The SO(4iV) GN model is a massive integrable field theory 
whose mass spectrum is known exactly .i^^jllij^ consists of the 
elementary fermions with mass m, their bound states, and of 
kinks. The bound states have masses {N > 1): 



m 



m 



^^'^ \2(2N-1)) 



with n 



(2(2JV-l)j 

2N — 2, while the kinks' mass reads: 



^^kinks 



2 sin 



i^ \2(2N-1)) 



(33) 



(34) 



The N = 2 case is special since the kinks mass is equal to 
that of the fermions. The SO(8) GN model enjoys a triality 
symmetry which has been exp loited in the study of the half- 
filled two-leg Hubbard ladder.l^^lSl In the iV > 2 case, the 
lowest excitations are the fermions which transform into the 
vectorial representation of the SO(4Af) group: they have the 
same quantum numbers as the original fermions R^ , L^ and 

The kinks transform in the spinorial representations of 
SO(4iV). It is more transparent to characterize these states by 
giving their charge and spin quantum numbers under U(l)c 
and SU(2A^) respectively: the kinks are those 2^^ states that 
carry charge Qk — k — N and transform in the ujk represen- 
tation of S\J(2N) where k varies from to 2N. One can dis- 
tinguish even and odd kinks, that transform in the even (odd 
respectively) spinorial representation and corresponds to even 
(odd respectively) fc's. In particular, the low-energy spectrum 
of the SO(4iV) GN model contains Sp(27V) spin-singlet states 
with charge Q — ±N which can be viewed as the generaliza- 
tion of the Cooperon excitations of the N — 2 case !— - I-^ These 



kink states identify with those discussed in Sec. II A 



The development of the strong-coupling regime in the 
SO(4iV) GN model leads to the generation of a spectral gap 
and the formation of a SP phase for all iV > 1 with order- 
parameter: 



O: 



SP 



(4i?„ - H.c 



(35) 



This order parameter is the continuum limit of the Spin Peierls 
operator Ogp = (—1)' 'J2a ^a i+i'^a.i and it has a non-zero 



expectation value in the GS as can be seen by a direct semi- 
classical approach of the SO(4:N) GN model. The phase is 
two-fold degenerate and breaks spontaneously the one-step 
translation symmetry (TaQ): La — > —iLa,Ra -^ iRa since 
Osp -^ -Csp under Ta„. This SP phase contains the \J{2N) 
line {V = 0) with C/ > 0, i.e. the repulsive \}{2N) Hubbard 
model. 

The second region of Fig.[3]can be easily determined with 
help of the duality symmetry Vli. The transformation of the 
SO(4iV) isotropic line under i^i is 1 = .92 /.9i = 93/91 = 
^94/91 = ~.95/ffi' which turns out to be the asymptote of 
the RG flow in the region (II) of Fig. [3] We thus deduce a 
second Mott-insulating phase which is obtained from the SP 
phase by applying the duality symmetry fii. Since La — > iLa 



under ili, its order parameter can be obtained from Eq. ( 35 1 



O. 



CDW 



LiRa+H.c. 



(36) 



which is nothing but the continuum limit of the 2kp CDW 
operator: Ocdw = Y.ai^^y^'L.i^aA- Region (II) is thus 
a fully gapped CDW phase. This CDW phase contains the 
U(2iV) fine (V = 0) with U < 0. 

What happens in region (III) of Fig. [3]is clearly of a differ- 
ent nature: the RG flow displays no symmetry enlargement, 
and we will have to develop other tools to tackle the physics 
in this interesting region. This will be done in the next section, 
where we will reveal striking differences according to the par- 
ity of N. Before that, we would like to give hand-waving 
arguments, based on the spectrum of the SO(4A^) GN model, 
supporting this even-odd scenario. To understand what hap- 
pens to the system when one leaves regions (I) and (II) where 
symmetry enlargement occurs, one should recast the whole 
particle content of the SO(4A^) GN model in multiplets of the 
internal continuous symmetry group of our problem, namely 
Sp(2A^)xU(l)c. One already knows how the SO(4iV) multi- 
plets split into U(l)c X SU(2Af) representations, that we write 
{Q, A) with A an SU(2iV) weight and Q the U(l)c charge (the 
number of fermions measured with respect to the GS). Denot- 
ing the vectorial representation (to which the GN "fundamen- 
tal fermions" belongs) by V, the even spinorial (to which even 
kinks belong) by 5'+-', and the odd spinorial (to which odd 
kinks belong) by 5*^^^ one has: 



V 

Si+) 



(l,Wi)® (-l,LJ27V-l) 

®f=o(2fc-7V,c^2fe) 
®^-Q{2k + l^N,uj2k+i)- 



(37) 



A quick way to check those quantum numbers is to note that 
they must be compatible with the fundamental fermions be- 
ing a boundstate of two kinks. The only missing piece of in- 
formation is the splitting of the SU(2iV) representations ujk 
into Sp(2iV) representation. Denoting by uik the fc* Sp(2A^) 
fundamental representation, one has: Ld2n = ®k=a^2k and 
'^2n+i — ffifc=o'^2fc+ii so that any SU(2A^) representation uj2n 
contains one and only one Sp(2Af) singlet, whereas all states 
of a;2„+i carry non-zero Sp(2Af) spin. 

Now one notices that region (III) of Fig.|3]occurs at negative 
V, where the system tends to favor Sp{2N) singlets. Let us as- 
sume that there is adiabatic continuity in the low-energy part 



of the spectrum. Then, the quantum numbers of the lowest en- 
ergy modes can be obtained by looking at those states in the 
SO(4iV) GN spectrum that are Sp(2N) singlets. It results that 
when N is even, one expects the "elementary" charged parti- 
cle (with the smallest U(l)c charge) to carry charge Q = ±2. 
On the other hand, when N is odd, there is a kink state that 
is a Sp(27V) singlet and carries charge Q = ±1. We will 
shortly see that this even-odd dichotomy does indeed occur, 
and that the elementary charged particles have the aforemen- 
tioned charges. 



D. Even-odd scenario 

The last region of the RG flow of Fig. [3] i.e. region (III), 
is difficult to analyze due to the absence of the second duality 
symmetry ^2 when N > 2. In this region, which includes 
the SU(2)c line with V = NU and U < 0, the operator with 
coupling constant gi in the low-energy effective Hamiltonian 
(24 1 reaches the strong-coupling regime before the others. In 



the limit of weak coupling, one has a separation of energy 
scales, due to the property ( [32] i of the RG flow. Neglecting 
all other couplings for the moment, the corresponding pertu r- 
bation is an integrable massive field theory for gi > QJ ^^ I ^^ lA 
spin gap As thus opens for the Sp(2N) spin sector in region 
(III). The next step of the approach is to integrate out these 
spin degrees of freedom to derive an effective Hamiltonian in 
the low-energy limit E -^ As from which the physical prop- 
erties of region (III) will be deduced. 



1. Parafermionization 

The resulting low-energy effective Hamiltonian involves 
the remaining degrees of freedom of the initial conformal 
embedding (19 1, i.e., the SU(2)7v sector Since the global 
continuous symmetry of model (24i is, in general, U(l)c x 



Sp(2A^), we need to understand how we go from the SU(2)Ar 
CFT to the U(l)c one. Such a mapping is realized by the 
conformal embedding: Zn ^ SU(2)Ar / U(l)c, which de- 
fines the Zat parafermio nic C FT series with central charge 
c = 2{N - 1)/{N + 2).EIlIlThis CFT describes the criti- 
cal properties of two-dimensional Zjv generalizations of the 
Ising model. The Zjv CFT is generated by the parafermionic 
currents 'ikL,R. {k — 1,...,N) with scaling dimensions 
hk = k{N~'k)/N. 

The different operators of Eq. ( [24| can be written in this 
parafermionic basis. First of all, the SU(2)7v currents (J2T|i can 
be directly expressed in terms of 'i'lL.R and a bosonic field 
$c which accounts for charge fluctuationsPl 



Jl,r - ^^«-P 



(±iy^87r/N 



$, 



cLM 



*1 



LM 



Jl 



L,R 



a.$ 



cL,R, 



(38) 



where the charge field <I>c = ^cL + ^cR is a compactified 
bosonic field with radius Re = yjNl2-n: <^c^ ^c^ \J2-eN . 



The remaining currents of Eq. ( 24 1 can also be expressed in 
terms of the parafermionic degrees of freedom using the re- 
sults of Ref. 23: 



^^ p 



/X2Tr(/)(2) exp UyJST:/N $c) 



(39) 



where 0^^' is the second primary operator of the Sp(27V)i 
CFT with scahng dimension 2N/{N + 2). In Eq. (|39]l, ei 
is the first thermal operator of the Z^r CFT with scaling di- 
mension A/{N + 2), and /i2 is the second disorder operator 
with scaling dimension 2{N — 2)/N{N + 2) which orders 
when the Zn symmetry is not spontaneously broken. ^^ 

Before investigating the low-energy limit E <C A^, it is 
crucial to analyze the hidden discrete symmetries of model 
(|24]| which become explicit thanks to the conformal embed- 
ding. It is well known that the Zat CFT has a global Zn x Zn 
discrete symmetry under which the parafermionic currents 
"ifkL (respectively '^kR) carry a (fc, k) (respectively (fc, —k)) 
chargep2] 

■^kL,R ^ e^^'^^'/'^^feL.ii under Z^, (40) 

with m = 0,...,A^— 1. The thermal operator ei transforms 
as a singlet under these discrete symmetries while the order 
and disorder operators ak, Hk carry respectively a (fc,0) and 
(0, k) charge: 



ak -> e*2^"'=/^crfc under Zn 



(41) 



and (Tfe (respectively ^k) remains unchanged under the Z^v 
(respectively Zm) symmetry. The Zm symmetry of the 
parafermions has a simple interpretation in terms of the orig- 
inal lattice fermions or the Dirac fermions of the contin- 
uum hmit. It is nothing but a special phase transformation 
or, in the continuum description: 






ia ^ e-'^^^/^L^, R^ -> e-''^'"/^i?„ 



(42) 



with TO = 0, . . . , A^ — 1. This Z^ symmetry leaves invariant 



model (24 1, and the correspondences (38 39 1 are also com- 
patible with the definition (42 1. In contrast, the Zat symmetry 
of the parafermions does not exist on the lattice. Away from 
half-filling, it becomes an independent emergent sy mmetr y of 
the model in the continuum limit and takes the formP323 



in 



e-""/^L„, Ro, -^ e"™/^i?„. 



(43) 



At half-filling, this transformation is no longer a symmetry of 
model (24 1 due to the umklapp operators. The Zjv symmetry 
has a more subtle role here: its combination with the following 
identification on the charge bosonic field: 



^r 



^r 



m 



INt: 

P\l ^r^ TO = 0, 



,7V-1, (44) 
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becomes a symmetry of model ( 24 1, as it can be seen from 
Eq. ( [39| ). In fact, this symmetry is a gauge redundancy since 
it corresponds to the identity in terms of the Dirac fermions. 
The last important discrete symmetries of the problem are the 



fli duality transformation ( 3 1 1 and the one-step translation in- 
variance (Tag), which only affect the charge field: 



r^i 


1 Nt' 


r 





(45) 



2. Low-energy Hamiltonian 

We are now in position to derive the low-energy limit E <^ 
As by integrating out the gapful Sp(2N) degrees of freedom. 



Using the parafermionization formulae (38 39 1, one finds 



'lint — 



+ ^(Ai2exp(iv/8^$,)+H.c.) 



(46) 



H.c, 



4,5 



(Tr0(2 



)52,4,5, and A3 = g3,N/2Ti. The 



where A 2 

low-energy Hamiltonian ( |46| ) enables us to explore the whole 
phase diagram of the model for all N . Along the SU(2)c line 
with V = NU, model ( |46l ) can be written in terms of the 
SU(2)Ar fields: 



H„„=53Ji?- JL + A2Tr$(i), 



(47) 



where $(^) is the spin-1 primary field of the SU(2)Ar CFT 
with scaling dimension 4/(iV + 2). The effective Hamil- 
tonian (47 1 is the low-energy theory of the spin-A^/2 SU(2) 
Heisenberg chain derived by Affleck and Haldane in Ref. [59] 



As shown by these authors, model (47i has a spectral gap, 
when N is even, while it describes a massless flow to the 
SU(2)i CFT when N is odd, in full agreement with Haldane's 
conjecture. ^^ The latter result has also be found by means of a 



parafermionic approach similar to Eq. (46 1 in presence of an 
SU(2) symmetry.''*^ 

The crucial point to map out the general phase diagram of 



the low-energy Hamiltonian ( 46 1 for all N stems from the sta- 
tus of the Zjv symmetry (42 1. The first term in Eq. (46 1 de- 
scribes an integrable deformation of the Z^v CFT which is 
always a massive field theory for all sign of A2.'^In our con- 
ventions, if A2 > (respectively A2 < 0) the Zat symmetry 
is unbroken (respectively spontaneously broken) and the dis- 
order fields (respectively order fields) condense: (/ifc) ^ 
(respectively {at) 7^ 0) for all fc = 1, . . . , iV. 

Let us first re-investigate the emergence of the CDW, SP 
phases in regions (I, II) within this parafermionization ap- 
proach. When A2 > (i.e. f/2 > 0), the Zat symme- 
try remains unbroken and one may integrate out the gapful 
parafermionic degrees of freedom to derive an effective field 
theory on the charge bosonic field. Since we have {112) 7^ 



and {^iL^iji! 7^ in the Z^v high-temperature phase, we 
obtain from Eq. ([46|: 



He 



Vc 

2 



Kr 



9c COS 



[V^tt/N 



$, 



where the Luttinger parameter is given by: 



Kr 



1 



^l + g3N/{2TrvF)' 



(48) 



(49) 



The low-energy Hamiltonian for the charge degrees of free- 
dom (48 I is the well-known sine-Gordon model at /3^ = 
SttKc/N. We thus deduce the existence of a charge gap when 
Kc < N which is always the case at weak coupling as seen 
from Eq. (49 1. The nature of the Mott-insulating phase de- 



pends on the sign of gc which is changed by the duality trans- 
formation fli (45 I. When gc < (i.e. 54 < 0), the develop- 



ment of the strong-coupling regime of the sine-Gordon model 
( |48| ) is accompanied by the pinning of the charged field on the 
minima: ($c) = 'P\J~Nt^J2, p being an integer. Since we have 
the identification ^c ^ ^c + V^ttN due to the periodicity of 
the charge field, we deduce that the insulating phase is two- 
fold degenerate with minima: (<I>c) = and ($c) = \/Nt:/2, 
i.e. the one-step translation symmetry (45i is spontaneously 
broken. The low-lying excitations are massive kinks and an- 
tikinks which interpolate between the two GS. The charges 
associated to these excitations are 



Q = ±^/2N/tt I dx d^^c = ±N. 



(50) 



For A^ = 2, the excitations correspond to the Cooperon 
excitations of the half-filled two-leg Hubbard ladder.!^ The 
charge excitations ( [SO] ) correspond to the generalization of 
these Cooperons. That they are the charge excitations with 
the minimal charge can be deduced from considerations on 
symmetry: amongst the spectrum of the SiO{AN) GN model, 
they are the only charged states that are both Sp(2iV) singlets 
and neutral under Zjv (indeed, the Zat charge of any state can 
be simply read off from the way it transforms under SU(2iV): 
states in w/j carry a Z^r charge given by k mod N). 

The physical nature of the two-fold degenerate Mott- 
insulating phases can be determined by expressing the SP and 
CDW ord er p arameters (35 36 1 in terms of the charge and the 
In fieldsP 

LlRc, ^ exp (i^2TT/N $,) A'lTr (j)'^^\ (51) 

where (/)(^) is the first Sp(2Af)i primary field with scaling di- 
mension (2iV + l)/2(iV + 2). Averaging over the Sp(2iV) and 
Ztv degrees of freedom, we obtain: 



CcDW 
Osp 



30S (V2vr/7V $,) 
3in (^/2tiIN < 



(52) 
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The phase with gc < (i.e. 54 < 0) is thus a CDW phase 
((CcDw) 7^ 0) and corresponds to the region (II) of Fig. [3] 
The second phase with g^ > is obtained from the CDW 
phase by the application of the duality transformation J7i. The 
pinnings of the charge field are then: (<l>c) = ^/NttJS and 
(<l>c) = 3-\/A^7r/8, which signals the formation of a SP phase 
in the region (I) of Fig. [3] since from Eq. (J52|, (Osp) 7^ 0. 
The quantum phase transition between the CDW/SP phases 
belongs to the self-dual manifold of the duality symmetry 51 1. 
Using the definition ( [45] l, one finds that the low-energy Hamil- 
tonian of the transition is given by 



n. 



^ If dual 



2 \Kr ' 



Kc {d^e. 



gc cos 



(V32^/iV 



$, 



(53) 



The resulting quantum phase transition is of BKT type. The 
transition displays a quantum-critical behavior with one gap- 
less bosonic mode if Kc > N/A. At this point, we need com- 
plementary numerical techniques to extract the value of Kc in 
order to conclude on the nature of the transition. 

Finally, the case with A2 < (or (72 < 0) corresponds to 
region (III) of Fig. [5] where the Z^r symmetry is now sponta- 
neously broken. In this Z^r low-temperature phase, the Z^r 
degrees of freedom are still fully gapped and the disorder op- 
erators now average to zero: (/ifc) = 0. Similarly to the 
A2 > case, we can integrate out the parafermionic fields 
to obtain an effective field theory on the charge bosonic field. 



However, due to the presence of the /X2 operator in Eq. (46 1, 
the resulting integration strongly depends on the parity of N . 
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FIG. 4: Values of the coupling constants gi in the infrared limit, in the 
vicinity of the transition from CDW phase (II) to BCS critical phase 
(III), close to the SU(2)c symmetric line (V^ = NU). Notations are 
the same as in Fig. [3] 



3. Phase diagram in the N odd case (N > 1) 

Let us first consider the case where N is odd. Since all 
the parafermionic operators in Eq. ( |46] l average to zero in 
the Zjv broken phase, one has to consider higher orders in 
perturbation theory to derive an effective theory for the charge 
field. The Zjy fields of model (46 1 carry a charge 2 under the 
Zjv symmetry (see Eq. (41 1 with k — 2 for /i2). When N 



is odd, one has to use the Nth order of perturbation theory to 
cancel out the Z^v charge of fj.2 so that we find: 



■H 



odd 



Vc 

2 



-^ (5,$c)' + Kc {d^Qcf 
Kc 



gc COS 



(VSnN 



$, 



(54) 



with gc ~ .94 , while we do not have any estimate of the 
Luttinger parameter except the bare one (49 1. On symmetry 



grounds, the effective Hamiltonian (54i can also be derived 
by finding the vertex operator in the charge sector with the 
smallest scaling dimension which is compatible with transla- 



tional invariance (45 1 and the redundancy (44 1. The resulting 



low-energy Hamiltonian ( p4| ) takes the form of a sine-Gordon 
model at /3^ = SttNKc so that a charge gap opens when 
Kc < 1/N. For Kc — l/N, this sine-Gordon model dis- 
plays a hidden SU(2) symmetry which should correspond to 



the SU(2)c line V = NU with U < that belongs to region 
(III). Close to this SU(2)c line, the RG flow of Fig. |4] shows 
that the coupling constant 94 is negative so that gc < 0. When 
Kc < 1/N, the charge bosonic field is thus pinned on the 
minima: ($c) = p-\/7r/2iV, p being an integer. Taking into 
account the gauge redundancy ( [44| , we find that the strong- 
coupling phase of the sine-Gordon model (|54]) is two-fold de- 
generate with (<!>£,) = and ($c) = ^7r/2A^. The charges of 
the massive kinks and antikinks excitations are now: 



Q = ±^/2N/Tr / dx d^^c = ±1, 



(55) 



in sharp contrast to the charge Q = ± A^ of excitations ( 50 1 of 



the CDW phase of region (II). At this point, we need to find a 
local order parameter to fully characterize the two-fold degen- 
erate Mott insulating phase in region (III). When the Sp(2A^) 
and Zat degrees of freedom are integrated out, the expression 



of the bilinear Dirac fermions ( 5 1 1 is naively short-ranged 
in region (III), since it contains the first disorder parameter. 



However, by fusing this operator with the Hamiltonian ( 46 1 
at the {N — l)/2 th order of perturbation theory, the disorder 
operator cancels out and one obtains the following low-energy 
description: 



LaRoc ~ exp 



(iV2nN 



$, 



(56) 
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FIG. 5: Values of the coupling constants gi in the IR limit, in the 
vicinity of the transition between the BCS (III) and SP (I) phases, 

displaying the restoration of an SU(2) symmetry, denoted SU(2)^, 
which is the dual of the lattice SU(2)c symmetry (V — NU): (72 = 
—(74 and g3 = —gs; we used the same notations as in Fig.lS] 



In region (III), where the Zjv is spontaneously broken, the 
CDW and SP operators then read as follows: 



CcDW ^ COS (V2'kN $c 
C^sp ^ sin (V2TrN $c) , 



(57) 



so that (C?cDw) 7^ 0. The insulating phase in region (III) 
when Kc < l/N is thus the continuation of the CDW phase 
of region (II). However, there is a striking difference at the 
level of the low-lying excitations: the generalization of the 
Cooperon excitations with charge Q = A^ is no longer a stable 
excitation in region (III) but becomes a diffusive state made of 



the kinks ( 55 1 which are massive holons. The situation is very 



similar to the SP phase of the half-filled U(4) Hubbard model 
between the weak and strong coupling regimes. ■'^ As far as 
the GS properties are concerned, there is a smooth crossover 
when the Z^r symmetry changes its status at A2 = and not a 
Zat quantum phase transition as it is the case away from half 

fiuingpmi 

When Kc > l/N, the charge degrees of freedom become 
gapless. We then need to determine the leading instability of 
this phase, i.e. the one that has the slowest decaying correla- 
tion functions. The singlet-pairing operator can be expressed 



in terms of the charge and the Zjv fields aspS 

Pjo - exp (iy/27:/N 6^) cnTr (^^^\ 



(58) 



Since the Zjy symmetry is broken, we have (<ti) 7^ and 
the low-energy representation of the singlet-pairing operator 

is thus: pJq ^ exp ( iy/2n/N 9c 1 . The gapless phase stems 

from the competition of this singlet-pairing operator, which 
cannot condense, and the CDW operator (57 1. The leading 



asymptotics of their equal-time correlation functions can then 
be straightforwardly determined: 



(Pjo(^)^oo(0))^Ax-V^^^= 



{n (x) n (0)) ^ -^ + i~ir^''° B x-^^^ (59) 

where n{x) is the continuum limit of the lattice density op- 
erator rii, and A,B are non-universal amplitudes. Since 
Kc > l/N, the leading instability of this gapless phase is 
the BCS singlet-pairing. 

The quantum phase transition between the gapful CDW 
phase and the gapless BCS phase occurs at Kc = l/N which 
coiTesponds to the SU(2)c (V = NU < 0) line. On this line, 
we observe that the exponents of the coiTelation functions of 
Eq. (59 1 are identical. Using the pseudo-spin operator ( 11 1, 



we deduce the following leading asymptotics from Eq. ( 59 1 



{sHx)s-m 
{s^x)s^m 



(-1) 



x/ao 



(60) 



The model with Kc = l/N displays a quantum critical behav- 
ior with central charge c = 1 for all odd TV and corresponds to 
the SU(2)i universality class. This result is in perfect agree- 
ment with the strong-coupling analysis of Sec. II along the 
SU(2)c line, where the pseudo-spin operator ( 1 1 1 is a spin- 



N/2, i.e. half-integer, operator. The low-energy properties of 
SU(2) half-integer Heisenberg spin chains are indeed known 
to be governed by the SU(2)i CFT.^'* In the spin language, the 
CDW and BCS phases are respectively the analog of the Ising 
and XY phases and the quantum phase transition occurs at the 
SU(2) Heisenberg point. 

When we deviate from the SU(2)c line in the gapless BCS 
phase. Figs. (|4]l5]) show that 53 decreases and then increases 
as a function of the interaction. Using the naive estimate of the 



Luttinger parameter (49 1, we deduce that Kc increases from 
Kc — l/N at the SU(2)c line and then decreases until one 
reaches the SP phase of region (I). The resulting transition and 
its properties can be deduced from the CDW/BCS transition 
by the duality symmetry J7i. Indeed, under the transformation 
( 45 1, the sign of gc of model (|54]) is changed and the gapful 
insulating phase when Kc < l/N is two-fold degenerate with 
($c> = V7r/8iV and ($c) = S^yn/SN. The SP order param- 
eter ( [57) 1 acquires a non-zero expectation value in this phase: 
(Csp) 7^ 0. As far as the GS properties are concerned, this 
phase is the continuation of the SP phase of region (I). The 
quantum phase transition between the BCS and SP phases oc- 
curs at Kc = l/N. Its position corresponds to an SU(2) line 
with g2 = —54, and 53 = —55 which is obtained from the 
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FIG. 6: Phase diagram obtained by the low-energy approach in the 
N odd case (A^ > 1). 



lattice SU(2)c line V = NU < (.92 = 54, 53 = 55) by the 
application of the duality symmetry fii. The resulting SU(2) 

line, noted 80(2)^ in Fig. pi does not exist on the lattice: it is 
an emergent SU(2) symmetry of the continuum limit. 

As a summary, Fig.l6]shows the zero-temperature phase di- 
agram of model (J2]i in terms of the lattice parameters U, V in 
the N odd case {N > 1), which results from the low-energy 
approach. 



4. Phase diagram in the N even case (N > 2) 

As in the N odd case, one has to consider higher orders in 
perturbation theory to derive an effective theory for the charge 



field $c since all the parafermionic operators in Eq. (46 1 av- 
erage to zero in the Z^v broken phase. When N is even, one 
needs the N/2 th order of perturbation theory to cancel out 
the Zn charge of fi2. The resulting low-energy Hamiltonian 
then reads as follows: 



n": 



2 
gc cos 






(V2nN 



$, 



(61) 



Alternatively, the effective Hamiltonian (61 1 can also be ob 



tained by considering the vertex operator in the charge sector 
with the smallest scaling dimension which is compatible with 



translational invariance ( 45 1 and the gauge redundancy ( 44 1 



The resulting low-energy Hamiltonian (61 1 takes the form of 
a sine-Gordon model at (3'^ — 2ttNKc so that a charge gap 
opens when K^. < A/N. One checks that, right on the SU(2)c 
line where the Luttinger exponent is constrained (Kc — 1/N), 
the sine-Gordon parameter takes the special value ^/2tt, at 
which it is known that a hidden SU(2) symmetry emerges.!^ 
The lowest energy modes are a massive triplet, the magnon 
of the integer spin Heisenberg model. Turning back to the 
generic situation Kc < A/N where a charge gap opens, the 
charge bosonic field is pinned into the following configura- 



if 5c < 




p being an integer The lowest massive excitations are the 
soliton and antisoliton of the sine-Gordon model; they carry 
charge 



g = ±^/2NJ^ I dx d^<^c = ±2, 



(63) 



which correspond to the Cooperon excitations. Using the 
gauge redundancy ( [44) i, we find that, in sharp contrast to the 
CDW and SP phases, the insulating phase when Kc < 4/7V is 
non-degenerate, its GS being: 



($,) =0, if 5c < 



($. 



, if 5c > 0. 



(64) 



Starting from the CDW phase of region (II), where the Z^ 
symmetry is unbroken, there is necessarily a quantum phase 
transition to the non-degenerate Mott-insulating phase of re- 
gion (III) with broken Zn symmetry. In particular, the dis- 
order parameter /ii of Eq. (J5T]| cannot be compensated using 
higher orders of perturbation theory as was the case for odd 
N. It means that in region (III): (Ocdw) ~ (C'sp) = when 
the Sp(2iV) and Z^ degrees freedom are integrated out. It 
is natural to expect that the non-degenerate insulating phases, 
described by the pinning ( [64| , signal the emergence of the HI 
and RS phases that we have identified in the strong-coupling 
approach (15i. At this point, it is worth observing that the 
duality symmetry ili plays a subtle role in the even N case. 



Indeed, under the transformation (45 1, the cosine term of Eq. 
( |6T] ) transforms as 



(65) 



cos (V27rN $c] -^ (-1)^/^ cos (V27rN $, 



so that there is room for an interesting N/2 even-odd effect. 

N/2 even case. 

Let us first consider the N/2 even case. A naive estimate of 
the coupling constant gc in higher orders of perturbation the- 
ory gives: gc ^ —54 . The RG flow close to the SU(2)c 
line {V = NU < 0) in the N/2 even case is similar to 
the one in Fig. |4] In this region, we have 34 < so that 
the non-degenerate gapful phase is described by the locking 
(<I>c) = of Eq. (J64|. As seen in Fig. [4] this region con- 
tains the SU(2)c line where the strong-coupling analysis ( 15 1 
predicts the emergence of the spin-iV/2, i.e. even spin, SU(2) 
Heisenberg chain. The low-lying excitation of the resulting 
Haldane phase is a gapped triplet state. From the expression 



of the pseudo-spin operator (111, one observes that it corre- 
sponds to a Cooperon excitation in full agreement with the 
prediction (63 1. We thus conclude that the Mott-insulating 
phase in the vicinity of the SU(2)c line (V = NU < 0), which 
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is described by the sine-Gordon model (61 1 with (<I>c) = 0, is 
the HI phase. 

The topological order of the Haldane phase with integer 
spin S = N/2 > 1 has been less understood than the S ~ I 
case. This phase displays edge states with localized spin N/A 
when OBC are used.^'' Unfortunately, we are not able to de- 
scribe these boundary edge excitations by means of our low- 
energy approach except when N = 2.'^^^**^^ On top of these 
end-chain states, the higher integer-spin Haldane phase should 
exhibit a non-local string ordering. ^^'™' A very naive guess is 
to use the generalization of the string-order parameter (28 i 
with spin-A^/2 operator In the low-energy limit, we find for 
N/2 even: 

lim (5fe"^^'-+i'5^'5|^) ~ 
K-i|-»-oo 



lim 

|2;-a|->-oo 



(sin ('v'^7r/2 $c (x)) si 
= 0, 



sm 



7V7r/2$e(y))> 

(66) 



since the HI phase is described by the pinning ($c) — 0. 
This result is in full agreement with what is known at the Af- 
fleck, Kennedy, Lieb, Tasak i (AKLT) point^' of the integer- 
spin Heisenberg chain,'^!!! and also from DMRG studies of 
the spin-2 Heisenberg chainl^S^A simple non-zero string or- 
der parameter in the HI phase, that we can estimate within our 
low-energy approach, is 



lim (cos 

\i-j\^OD 



V k<i / \ k<j 




lim (cos 

\x-y\-^oci 



i (V^7r/2 $e (x)) cos (^^NTr/2 $, (t/))) 
^ 0. (67) 



This lattice order parameter turns out to be non-zero at the 
AKLT point of even-spin Heisenberg chains. ^^ 

In summary, when N/2 is even, the HI phase is described 
at low-energy by the sine-Gordon model (61 1 with Kc < ^/N 
and a non-degenerate GS ($c) = 0. The quantum phase tran- 
sition between the CDW and HI phases is difficult to deter- 
mine exactly. On general grounds, we expect an Ising quan- 
tum phase transition or a first-order one due to the difference 
of the GS degeneracies between the two phases. In the CDW 
and HI phases, the charge bosonic field is locked at ($c) = 
so that the CDW/HI quantum phase transition is governed by 
the Zjv interacting field theory: 



nt 



X2 ei + A4 (/i2 +H.C.) . 



(68) 



Model (roSjl is a deformation of the Z^r CFT perturbed by two 
relevant operators with scaling dimensions 4:/{N + 2) and 
2{N — 2)/N{N + 2) respectively. When acting separately, 
each perturbation yields a massive field theory, but the inter- 
play between them may give rise to a second-order phase tran- 
sition at intermediate coupling. In this respect, when A2 < 0, 
the first operator in Eq. ( 68 1 orders the Z^ degrees of free- 
dom while the second one wants to disorder them. We con- 
jecture that this competition for A2 < leads to a massless 



flow to a Z2 quantum critical point in the IR limit. The quan- 
tum phase transition between the CDW and HI phases thus 
belongs to the 2D Ising universality class with central charge 
c = 1/2. In the simplest iV = 4 case, we can show this re- 
sult explicitly by exploiting the fact that the Z4 parafermionic 
CFT has central charge c = 1 and so it should be possible 
to realize it with a single free Bose field. In fact, the correct 
identification is quite subtle and the Z4 CFT turns out to be 
equivalent to a Bose field living on the orbifold line at radius 
R = y^3/27r.'22i However, as shown in the Appendix of Ref. 
l74l it is still possible to bosonize some fields of the Z4 CFT 
with a simple (periodic) Bose field $ defined on the circle with 
radius R = ^3/27r: $ ~ $ + 2nR. In this respect, the two 
operators of Eq. (681, with scaling dimension 2/3 and 1/6, 
take the form of vertex operators. The bosonized description 
of the effective field theory ([68]l reads: 



nt,t = A2 cos 



(y/8TT/3 $") + A4 sin (y/2^ 



$ 



(69) 



This model is the so-called two-frequency sine-Gordon model 
which, for instance, governs the transition from a band insu- 
lator to a Mott insulator in the ID ionic Hubbard model.l^Sl 
When A2 < and for all signs of A4, model ( [69| displays a 
Z2 quantum critical point in the IR limit which has been ana- 
lyzed non-perturbatively in Refs lTsHTTl We thus deduce that 
the quantum phase transition between the CDW and HI phases 
for iV = 4 belongs to the 2D Ising universality class. 

Let us now investigate the fate of the HI phase as one de- 
viates from the SU(2)c line. As in the N odd case, there is a 
regime in region (III), away from the SU(2)c line, where the 
coupling 53 that appears in the Luttinger parameter expression 
( |49l ) decreases and then increases as function of the interac- 
tion (see Figs. [4J [5]). In the vicinity of the minimum of g^, 
we expect the emergence of a gapless phase associated to the 
sine-Gordon model (61 1 with Kc > A/N . The existence of 



this intermediate gapless phase will be confirmed numericaUy 
in Sec. VI by means of DMRG calculations. In this respect, 
the N ^ 2 case is very special since this phase shrinks to 
a line which marks the phase transition between HI and RS 
phases (see Fig. |2]). This critical phase has only one gapless 
charge mode and the singlet-pairing has the same low-energy 



behavior as in the N odd case: P, 



00 .^ exp {i^27r/N O^j 

However, this phase is different from the gapless BCS phase 
of the N odd case. Indeed, as already stressed, the disorder 



parameter iii of Eq. (51 1 cannot be compensated using higher 
orders of perturbation theory which means that the alternating 
part of the CDW operator is short-ranged. We then deduce the 
following leading asymptotics of the equal-time correlation 
functions: 



{n{x)n{0)) --^, 



(70) 



where in the density correlator, only the uniform part has a 
power-law decay. The leading instability is the singlet-pairing 
one when Kc > A/N . The main difference with the gap- 
less BCS phase of Fig. l6] stems from the fact that the alternat- 
ing part of the density correlator (70 1 has now an exponential 
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decay. The quantum phase transition between HI and BCS 
phases belongs to the BKT universality class. 

The last regime of region (III), which corresponds to the 
transition between (III) and (I) of Fig. E\ can be identified by 
means of the duality symmetry fii. Under this transformation 
(|45]l, the CDW phase is changed into the SP phase. In contrast, 



from Eq. (65 i, we deduce that the sine-Gordon operator of the 



low-energy Hamiltonian (61 1 remains invariant when N/2 is 

even. The Mott-insulating phase in the vicinity of the SU(2)^ 
line of Fig.|5]is thus described by the strong-coupling regime 
of the sine -Gordon model (61 1 with K^ < 4:/N and the pin- 
ning (^c) ~ 0. This phase is expected to be the RS phase, 
i.e. the so-called large D phase of the integer spin Heisen- 
berg chain,'*'^ which appears in the strong-coupling approach 



( 15 I for a sufficiently strong positive D. Interestingly enough. 



within our low-energy approach, this phase is described ex- 
actly in the same way as the HI phase. Thus, the two phases 
necessarily share the same order parameters, such as the string 



orders (66 67 i for instance. However, they should have dif- 
ferent edge states but we could not, very unfortunately, inves- 
tigate these boundary end excitations in our CFT approach. 
Recently, it has been argued that the edge-state structure of 
the even-spin Heisenberg chain is not protected by symme- 
try in contrast to the odd case .1^ In particular, the authors of 
Ref "SThave conjectured that there is an adiabatic continu- 
ity between the Haldane and large D phases in the even-spin 
case. The Haldane phase is thus equivalent to a topologically 
trivial insulating phase in this case. This adiabatic continuity 
has been shown numerically in the spin-2 XXZ Heisenberg 
chain with a single-ion anisotropy by finding a path where the 
two phases are connected without any phase transition.lSl In 
our problem, the HI and RS phases are separated by an inter- 
mediate gapless BCS phase. However, within our low-energy 
approach, the two non-degenerate Mott-insulating phases are 



described in the same manner by the sine-Gordon model (61 1 
with Kc < 4:/N and the pinning ($c) = 0. In this respect, 
our results strongly support the conjecture put forward in Ref. 
[33] 

Finally, the quantum phase transition between RS and SP 



phases is obtained from the effective theory (68 i by the ap- 
plication of the duality symmetry J7i. The latter transforma- 
tion changes the sign of the coupling constant A4 of the ^2 
operator. However, this sign is irrelevant for the competition 



between the two relevant operators in model (681. We thus ex- 
pect that the resulting quantum phase transition still belongs 
to the 2D Ising universality class. In summary, Fig.lTlpresents 
the phase diagram, in terms of the lattice parameters U, V in 
the N/2 even case. 

N/2 odd case. 

The last case to consider is the case where N/2 is odd. 
In region (III), in the vicinity of the SU(2)c line, the charge 
bosonic field $c of the sine-Gordon model ^6T\ with Kc < 
4:/N is now pinned at ($c) = \Jt^ /2N since g^ > 0. The 
non-degenerate Mott-insulating phase is the HI phase. This 
phase can be described by the generalization of the string- 
order parameter (28i with spin-A^/2 operator Indeed, in the 




FIG. 7: Phase diagram obtained by the low-energy approach in the 
TV even case (iV > 2). 



low-energy limit and for the N /2 odd case, we find: 



lim (5fe"^'^=-+i'5^'5" 

|i-j|^oo 



lim 

\x-y\-^oo 



(sin (^JN■K/2 $e (a:)) sin (^N^ $c (y))) 
y^O, (71) 



in sharp contrast to the result ( 66 1 of the N /2 even case. For 



general odd-spin Heisenberg chain, the order parameter (|7^ 
is known to be non zero contrarily to the even-spin case.^^ 
In this respect, there is thus a clear dichotomy in the HI phase, 
depending on the parity of N /2. For odd-spin Heisenberg 
chains, the authors of Ref. [33] have predicted that the Hal- 
dane phase displays a topological order and is not equivalent 
to the large D phase as in the even-spin case. This scenario is 
in perfect agreement with our low-energy approach. Indeed, 
according to Eq. ( [65| , the duality symmetry J7i changes the 
sign of the vertex operator of model (61 1 when N /2 is odd. 
The physical properties of the RS phase are thus governed by 
the sine-Gordon model ( [6T] ) with K^ < 4/iV and the pinning 
((j)^) = 0. In the N/2 odd case, the HI and RS phases are 
described by two different locking of the charge bosonic field 
in sharp contrast to the N/2 even case. In particular, the RS 



phase is described by the string-order parameter ( 67 1 and not 
( [7T| ) as the HI phase is. The HI and RS phases are thus totally 
distinct phases that cannot be adiabatically connected. 

Finally, as in the N/2 even case, the transition between 
these two non-degenerate phases is accompanied by the for- 
mation of an intermediate gapless BCS phase with the prop- 
erties (70 1. Unfortunately, in the N/2 odd case, we do not 
have access to a theory of the quantum phase transition be- 
tween CDW (respectively SP) phase and the HI (respectively 
RS) phase. We suspect, as in the N/2 even case, that the tran- 
sition belongs to the 2D Ising universality class but it certainly 
requires a proof. Fig.lTlpresents the phase diagram in the N/2 
odd case which, apart from the subtleties on the topological 
nature of the HI phases, is identical to the N/2 even case. 
Last, we would like to emphasize that the N = 2 case (see 
Fig.l2]) is not representative of the even family but turns out to 
be special. 
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IV. PHASE DIAGRAM OF HALF-FILLED SPIN-3/2 
FERMIONS (iV = 2) 

In this section, we give the phase diagram of model ([T]l 
when A^ = 2 in the {U /t, V/t) plane, obtained from numeri- 
cal calculations. Four phases are found and reported on Fig. [8] 
two phases which break translational invariance, the SP and 
CDW phases, and two with non-degenerate GS which can 
only be distinguished through non-local string orders, the HI 
and RS phases. These phases are separated by transition lines 
determined numerically (full lines), and compared to weak- 
and strong-coupling predictions displayed in dashed lines. In 
addition, three particular lines are plotted where the model has 
an exact enlarged symmetry that we have discussed in Sec. 11. 

The numerical calculations are performed with DMRG on 
chains, each site containing the 16 states of the onsite basis 
(for N > 2, since the local Hilbert space is too large, we 
must use other strategies as discussed below). We fix three 
quantum numbers: the spin part 5^ = | X]q i(~l)"^^'^a.i' 



^z 



iEi('^i, 



"-2.1 



na,, 



n4_i) as well as the total 



number of particles Nf = 2L, i.e. the total charge. The GS 
lies in the S^ — T'^ — sector. The number of kept states 
is typically m = 2000 and OBC are used if not stated other- 
wise. Denoting by L the length of the chain, the local order 
parameters are computed numerically by taking their value in 
the bulk of the chain (we choose to work with an even number 
of sites): 



CcDw(i) = "-L/2 - "L/2-1 
Osp(i) = iL/2-<L/2-l, 



(72) 
(73) 



where Uj = J^a "« j i^ '■^^ '^^'^^^ onsite density and tj - 
J2a ^a j+i'-a 3 + ^.c. the local kinctic energy on bond (j, j 



A. The HI phase 

We start a more detailed discussion of the phase diagram 
from the V — 2U line which shows the remarkable SU(2)c 
symmetry, leading to the effective spin-1 Heisenberg model 



( 15 I in terms of charge degrees of freedom. We have recently 
demonstrated'^ that the gapped HI phase is realized for a given 
value of V/t, and that its extension is rather small. We here 
refine the description of the boundaries of the Haldane phase 
and discuss the nature of the transition lines to respectively the 
CDW and RS phases. In order to find the transition line from 
CDW to HI, we use Ocdw which vanishes in the HI phase 
and which is straightforward to compute. The transition from 
HI and RS is more difficult to determine as no local order pa- 
rameter can discriminate between the two phases. In Ref. [15] 
we gave several signatures of the transition which can be used 
to locate it: non-local charge string order parameters and the 
presence of edges states, which are observed here by consid- 
ering a charge excited state with two additional fermions; for 
OBC, this state has a vanishing gap to the GS. The simplest 
way to determine the transition with our numerical scheme is 
to look at the distribution of the charge in the excited state 




FIG. 8: (Color online) Phase diagram of the N = 2 case. Full lines 
are DMRG results (see text for discussion) and dashed lines are re- 
sults from solving numerically the RG flow in the weak-coupling 
limit, c stands for the expected central charge on the transition 
line. Lines with higher symmetries are indicated in black dashed 
lines. Boundaries of the HI phase obtained by the strong-coupling 
are shown in dotted lines. As a rule of thumb, DMRG results can 
be trusted if |(7|, | V| > i, while weak-coupling predictions are exact 
close to the free fermions limit at the origin. Inset: Zoom on the re- 
gion of the HI phase where the x-axis {2U ~ V)/t is perpendicular 
totheSU(2)c xS0(5)line. 



with Nf ^ 2L + 2: an excess Nf — 1 charge will be stuck 
at each edge in the HI phase (equivalent to the spin- 1/2 edge 
state of the spin-1 Haldane phase), while an Nf = 2 excita- 
tion lies in the bulk of the RS phase (equivalent to the 5=1 
magnon of the Heisenberg ladder). We thus use this change 
in the density profile of the charge excited state (benchmarked 
with other signals of the transition for V = — 2t) to give the 
estimate of the transition line in Fig. [8] 

In the weak-coupling regime, |[/|, |T^| < t, DMRG calcu- 
lations become hard due to the relevance of many low -energy 
onsite states. In the strong-coupling limit (large U, V), on- 
site energy scales are well-separated so that DMRG efficiently 
eliminates high energy irrelevant states. The numerical pre- 
dictions of the RG flow provides a better prediction for the 
transition lines in this weak-coupling regime: these estimates 
are V ~ 3.33C/ for RS-HI and V = 1.56U for CDW-HI. 



The two transition lines in Fig.lS 



are also compared to the 
strong-coupling predictions of Sec.lIIB For large \U\/t and 
|y|/i, the effective Hamiltonian around the SU(2)c line is a 
spin-1 model with antiferromagnetic coupling J = 2t^/5|C/| 
and anisotropy D = 2U—V (see Eqs. ( 15 16 1). The ph ase di- 
agram of this model has been extensively studiecPE^EMlol and 
shows that a Haldane-Neel transition (equivalent to the HI- 
CDW one) occurs for D/J ~ —0.5 while a Haldane-large- 
D transition (equivalent to the HI-RS one) is obtained for 
D/J ~ 1. This gives the two curves Vhi-cdw = 2f/ + t'^ /\U\ 
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FIG. 9: (Color online) Behavior of the Haldane gap A along the 
SU(2)c line V = 2U < of Fig. [8] (a) gaps and finite-size scal- 
ings (see text for discussion), (b-c) the Haldane gap as a function 
of U/t in units of respectively t and the effective antiferromagnetic 
coupling J. The As=i line indicates the value of the gap known for a 
Heisenberg spin-one chain. 



and Vhi-rs = 2f/ — 2t^ /\U\ explaining both the shrinking and 
the asymmetry of the extension of the HI phase in the strong- 
coupling regime. 

Although the Haldane gap decreases in the strong-coupling 
regime simply because J decreases, the agreement between 
the fermionic spin-3/2 Hubbard model under study and the 
spin-1 effective model improves as irrelevant degrees of free- 
dom are pushed to high energies. This can be illustrated nu- 
merically by the behavior of the Haldane gap along the SU(2)c 
line as a function of U/t. The Haldane gap is computed using 
OBC from the following gaps: 



A„b = Eo{2L + 26) - Eo{2L + 2a) 



(74) 



where EQ{Nf) stands for the GS energy with Nf fermions. 
As evoked previously, the presence of edge states with OBC 
makes the first excited state collapse onto the GS, so that Aqi 
vanishes in the thermodynamical limit. Still, both A02 and 
A12 must remain finite and tend to the bulk Haldane gap for 
sufficiently large sizes. These behaviors, together with finite- 
size extrapolations of the gaps using the ansatz 



^ab{L) = Aab(oo) + const. 



-^/Vi 



(75) 



are clearly shown by the numerical results of Fig. l9la). 
Fig.|9|b) and (c) display the extrapolated gaps as a function of 
U/t in units of respectively t and J. While the weak-coupling 
opening of the gap cannot be reliably studied here, we observe 
that the gap passes through a maximum around U/t ~ — 1 
which is close to value U/t ~ —1.5 for which the width 
of the HI phase is maximal. In the strong-coupling regime, 
the gap in units of t decreases as expected (see Fig. l9|b)), 
while, put in units of J (see Fig. |9lc)), it eventually reaches 
the value As=i ~ 0.41J known for the spin-1 Heisenberg 
chain: U ~ —4t is already deep in the strong-coupling regime 
along this SU(2)c line. 



cq 



HI-CD W, U/t = -1.04 



HI-RS, U/t = -0.89 
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FIG. 10: (Color online) Fitting the entanglement entropy close to the 
critical lines surrounding the Haldane phase provides central charges 
c close to the expected values c = 1/2 (for HI-CDW) and c = 1 (for 
HI-RS). The best agreement is found using PBC with DMRG and 
keeping a large number of kept states m. 



Lastly, we investigate the nature of the critical points at the 
two boundary lines of the HI phase. From the low-energy re- 
sults of Sec. Ill B, we expect that the Hl-CDW is an Ising tran- 
sition with a central charge c = 1/2 while the HI-RS transi- 
tion belongs to the BKT type, associated with a central charge 
c = 1. In the strong-coupling limit, this has been observed 
numerically for the spin-1 chain with single-ion anisotropy.l^ 
To check these predictions from the DMRG data, we use the 
universal scaling of the entanglement entropy (BE) in a critical 
phase, which gives a direct access to the central charge. We 
obtained the most convincing results using periodic bound- 
ary conditions (PBC) at the price of keeping a much larger 
number of states and using small system sizes. Similar cal- 
culations have been peiformed in the context of the SU(A^) 
generalization of Haldane's conjecture-l^ The results of the 
EE on a finite chain of length L — 48 along the V = —2t 



around the HI phase are given in Fig. 10 The central charge 
is obtained from the data using the universal formula^ 



S{x) = - lnd{x\L 



1) + const. 



(76) 



with d{x\L) = 7 In (^) the cord function and S{x) the EE 
of a block of size x with the rest of the chain. The values 
obtained for c are in good agreement with the expected val- 
ues considering the large number of local degrees of freedom. 
There is an uncertainty on the location of the critical points 
but, on a finite system, as long as L <C ■C, with ^ the correla- 
tion length associated to the closing gap, the physics will be 
effectively that of the critical point. 



B. The RS-SP transition 

We now turn to the discussion of the RS-SP transition in 
the right-down quadrant of Fig. |8] The two phases RS and 
SP can be simply distinguished by the local spin-Peierls order 
parameter C^sp which is finite in SP while it is zero in RS. The 
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FIG. 11: (Color online) Scaling of the order parameter Ogp at the 
RS-SP transition : (a) for U /t = 4, the exponent is quite close to 
1/8; (b) for U/t = 1, it is closer to 0.8. 



vanishing of the order as the system size increases provides a 
good estimate of the transition Hne. 

We further try to give evidence for the nature of the transi- 
tion and check whether it lies in the Ising universality class. A 
possible approach is to use the EE again and look for c = 1/2. 
However, the Osp order parameter appears as the leading cor- 
rections to the EE with OBC and gives strong oscillations in 
the signals, particularly in the SP phase and up to the criti- 
cal point. These oscillations render the fits difficult, and the 
value of c is not reliably extracted for the accessible system 
sizes. Using PBC improves a bit the situation, but the oscillat- 
ing parts of the EE could not be suppressed (as expected for 
the GS) with DMRG, even by increasing the number of kept 
states and sweeps. 

Consequently, we use another strategy to identify the Ising 
universality class. We know that the correlation function of 
the order parameter has a universal exponent 1/4 at the crit- 
ical point. Then, Friedel oscillations of the order parameter 
gives the scaling C'gp(-L) ex L^^l^ on the critical point. In 
the SP phase, Osp(i) reaches a constant in the thermody- 
namic limit, while it decreases exponentially in the RS phase. 
Thus, by looking at the scaling of C'sp(_L) for different pa- 
rameters, we are able to give both a precise estimate of the 
critical point and to check that the exponent is indeed close 
to 1/8. The results along two cuts at [/ = 4t and [/ = i are 
reported in Fig. 
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In the strong-coupling regime U = 4t, 
we do observe a very good agreement with an exponent 1/8, 
typical of the Ising universality class. However, in the weak- 
coupling regime, a much larger exponent of 0.78 ~ 6/8 fits 
well the scaling curves. We understand this discrepancy in 
the following way: in the weak-coupling regime, the gaps to 
higher excited states are too small to be thrown away in the 
low-energy regime of ■a. finite system. In other words, the cor- 
relation lengths associated with theses gaps become too large 
and we could not reach sizes sufficiently large to freeze them. 
A speculative picture can account for the observed number: 
at weak coupling, the transition line gets very close to the 
SO(6) line which has the equivalent of six gapped Ising de- 
grees of freedom, but with an exponentially small gap of the 



order i/C/.I^In this weak-coupling regime, the numerics can- 
not resolve these gaps and the Ising degrees of freedom appear 
critical, each contributing to 1/8 in the exponent which then 
should be close to 6/8. 

This comment brings us to the discussion of effect of the 
proximity of the SO(6) line (an exact enlarged symmetry) to 
the RS-SP transition line. The y = and t/ > line has been 
studied analytically and numerically in Ref. 37 the charge 
and spin gaps open slowly with IJ jt and are numerically neg- 
ligible below I] ~ 2<. In the weak-coupling regime, the low- 
energy physics has an emerging enlarged SO(8) symmetry. In 
the strong coupling regime, the spin gap decreases after pass- 
ing through a maximum around U ~ 6f . The data shows that 
the RS-SP transition line has a non-monotonic behavior, first 
following the weak-coupling RG predictions and then being 
attracted by the SO(6) line at large interactions (see Fig. [Sll. 
This attraction can be qualitatively understood by the behav- 
ior of the spin gap as \J increases. Considering y as a pertur- 
bation which closes the spin gap A5, the line should typically 
behave as VAU) ^ — A5(L/) which is non-monotonous and 
stick to the SO(6) line in the strong-coupling limit. In the 
weak-coupling limit U < t, the RG prediction V = —O.IOU 
is more reliable than the numerics. 



C. The CDW-SP transition 

Lastly, we briefly discuss the CDW-SP transition between 
these two phases which breaks translational symmetry. Nu- 
merically, the precise determination of the transition with Ogp 
and OcDW using DMRG turns out to be difficult due to for- 
mation of domains of each kind of orders close to the transi- 
tion line. Changing the number of kept states, the number of 
sweeps and the size, slightly moves the transition point deter- 
mined by the order parameter at the center of the chain. This 
leads to error bars in the phase diagram which are relatively 



small compared to the parameter scales of Fig. 1 1 but are too 
large to focus on the critical features of the transition line. We 
could not check the c = 1 expectation of this transition, due to 
both the difficulty in locating the transition point, and because 
of strong SP oscillations in the EE. Notice that on the criti- 
cal line, the correlations of the quartet operator c| ^c^ jCg jcjj j 
become critical which is qualitatively in agreement with nu- 
merical observations. 

Here again, we see that the transition line is rather close to 
a high symmetry line of the phase diagram, namely the SO(7) 
line V = — 2C/.^' In the weak-coupling regime, the numerical 
solution of die RG Eqs. (l30| for iV = 2 gives V = -1.61U 



but, for larger |C/|, DMRG calculations indicate that the tran- 
sition is attracted to the vicinity of the SO(7) line. A argument 
similar to the one used for the RS-SP transition can be drawn: 
we see that SO(7) line is in a SP gapped phase. The strong- 
coupling spin-model along this line is an SO(7) Heisenberg 
model where the spins belong to the vectorial representation 
of SO(7)^' and our analysis predicts a SP bond ordering. Nu- 
merically, we compute the spin gap Ag and charge gap Ac 
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FIG. 12: (Color online) Charge and spin gaps and their ratio along 
the S0(7) line V = -2(7 of Fig.[8] 
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FIG. 13: Numerical phase diagram obtained by DMRG in the TV = 3 
case. 



defined as follow: 

A, = EoiN,l) + EoiN,-l) -2Eo{N,0) 

A, = Eq{N + 2,0) + EoiN -2,0) -2Eo{N,0) 

where Eo{Nf, S^) is the GS energy with Nf fermions in the 
S^ sector with T^ = and iV = 2i is the reference number of 
particle at half-filling. The results extrapolated in the thermo- 
dynamic limit are given in Fig. 12 for a wide range of U/t val- 



ues. The gaps open slowly in the weak-coupling regime and 
then reach a maximum around C/ ~ —7t, before decreasing in 
the strong-coupling regime. The ratio of the gaps A^/A^ is 
very close to two, everywhere but in the weak-coupling limit 
where the numerics are challenging for accurate predictions. 



V. PHASE DIAGRAM IN THE iV = 3 CASE 

In this section, we investigate the phase diagram of 
model ([T]| when A^ = 3 and in the {U/t, V/t) plane using ex- 
tensive DMRG simulations. Since the local Hilbert space on 
each site contains 2^ = 64 states and is quite large, we have 
implemented the following strategy: we use a mapping to a 3- 
leg Hubbard ladder where the chains correspond to fermionic 
states with Sz equal to ±1/2, ±3/2, and ±5/2 respectively. 
Then, after some algebra, we can rewrite all hoppings and in- 
teraction terms in this language, which introduce for instance 
rung interactions and rung pair- hopping terms. This mapping 
to a ladder allows us to converge faster to the GS, but we have 
checked that the symmetry between chains is preserved in the 
SU(3) case for instance. Typically, we keep between 1600 and 
2000 states in our simulations for measuring local quantities 
and up to 3000 for correlations, and we use OBC. 

Since no topological phase is expected, we can rely on mea- 
suring local quantities such as local density and kinetic energy, 
as well as density and pairing correlations that will character- 
ize the critical phase that has been shown to exist along the 
SU(2)c line V = W in Ref. JS, The following phase di- 



agram can thus be obtained in Fig. [13] and it contains only 
three phases: SP, gapless BCS and CDW. 

Data points on this plot correspond to simulations done on 
system length L — 72, while phase boundaries were also ob- 
tained from scaling different system sizes (see below). 



A. Properties along the SU(2)c line 

We start by considering the SU(2)c line V — 3U. For large 
enough \U\/t, the strong coupling argument of Sec. II B tells 
us that the chain will behave effectively as an antifeiTomag- 
netic Heisenberg spin-3/2 chain, which is known to be critical. 
In Fig. [I4{a), we show how pairing and density correlations 
behave along this SU(2)c line. Their long-distance form has 



been determined in Eq. ( 59 1 and, measured from the middle 
of the chain, reads: 



V{x) 
J^{x) 



{PUL/2 + x)PooiL/2)) 



A 



(77) 



{n{L/2 + x)n{L/2)) - {n{L/2 + x)){n{L/2)) 



NKc {-l)''B 






rNKa 



Using the definition of the pseudo-spin operator ( [TT| , we ob- 
serve that the two coiTelations match perfectly, as expected of 
course for an exact SU(2)c symmetry. Both correlations are 
algebraic and expected to decay as Vlna;/a;,l^but it is known 
that checking accurately logarithmic corrections is a challeng- 
ing numerical problem**** that we will not further investigate. 

Another peculiar property of spin-3/2 chain with OBC was 
conjectured by Ng,'^and confirmed later numerically Pleven 
though the system is critical, one can observe "edge states" 
with OBC, in the sense that the magnetization profile will ex- 
hibit an excess close to the edges, although there are no finite 
coiTelation length (i.e. the magnetization profile decays alge- 
braically away from the edges). Here, we investigate a similar 
situation, namely with a charge SU(2)c symmetry where it is 
the local density that plays the role of the magnetization for 
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FIG. 14: Pairing and density correlations obtained by DMRG in the 
A'^ = 3 case for various interactions corresponding (a) to the ex- 
act SU(2)c symmetry; (b) to the emergent SU(2) symmetry SU(2)^. 
Note that correlations are measured starting from the middle of the 
chain. 



actual spin-3/2 chain. When adding 2, 4 or 6 particles (with 
respect to half-filling), as shown in Fig. 



15 we do observe 



modulations in the local densities reflecting these edge states. 
Physically, it means that the first excitation (adding 2 parti- 
cles) is an edge excitation, but the next ones correspond to 
making a bulk excitation. 



iVf = 3L + 6 



{"(l)>jV/.3i-2 - 3 

("(')>N/=3L + 4 - ("W)K/«3i ■ 2 




FIG. 15: Local densities obtained by DMRG in the A'^ = 3 case 
for U/t = -2, V = 3(7, and L = 108. From bottom to top, 
data correspond to adding 2, 4 or 6 particles to the half-filled system. 
Data for adding 4 and 6 particles are shifted by 2 and 4 respectively 
for clarity, and in these cases {n{i))Mf=3L+2 has been substracted 
in order to get the bulk contribution. 



large extension. As will be shown below, this critical phase 
has dominant BCS pairing correlations, thus its name. For 
fixed negative V, we observe the transition to the SP phase 
for large enough U > 0. This is in agreement with the low- 
energy prediction and the RG phase diagram (see Fig. [6]). 

In order to characterize the critical phase, we can compute 
its Luttinger parameter Kc from the behavior of either pairing 
or density correlations, using Eq. dTTJ). In Fig. 16 we plot both 
correlations at 1/ = — 6i and for various values of U. Let us 
start with the discussion of Vix) (which coiTesponds up to a 
factor 2/3 to the transverse pseudo-spin correlation function). 
In order to be able to fit over the whole range, '^^ data are plot- 
ted vsx' = d{x\L+l) / y/cos{TTx / {L + 1)), where d{x\L+l) 



is the cord function, defined in Eq. ( |76] l. We observe a very 
smooth behavior, which allows to extract the behavior of Kc 
vs U (see Inset). Due to the logarithmic corrections which are 
known to exist along the SU(2)c line, it is very hard to recover 
that Kc —> 1/3 when U -^ — 2t as expected from the exact 
SU(2)c symmetry. Moving away from the SU(2)c line, our 
data indicate that Kc rapidly reaches a maximum, before go- 
ing down again. The transition to SP coiTesponds to pairing 
correlations that become exponential (not shown), and occurs 
when Kc = 1/3 in agreement with the low-energy approach. 



Another way to compute Kc consists in using Eq. ( 77 1 for 
density correlations. In principle, one can use either the uni- 
form or alternating part to extract it. However, in the regime 
where Kc < 2/3, the alternating pait is dominant, whereas 
in the opposite case, the uniform part decays more slowly. 
Therefore, we have fitted either the uniform part or the stag- 
gered part to extract the value of Kc shown in the Inset. 

Overall, we have an excellent agreement between the es- 
timates of Kc obtained from both correlations, which gives 
confidence in the validity of the Luttinger liquid description of 
this critical phase. Moreover, the behavior of Kc vs U is com- 
patible with our expectation (see Sec. |IIID3] i: Kc exceeds 1/3 
in the BCS phase (giving rise to dominant BCS correlations) 
and the transitions to SP and CDW occur when Kc reaches 
1/3. 




B. The transition from critical BCS to SP 



As can be seen from the phase diagram shown in Fig. 13 
the critical phase that exists along the SU(2)c line has a rather 



FIG. 16: (Color online) Pair and density correlations in the A'^ = 3 
case for V/t = — 6 and L — 72 and various U/t. Inset: Fitting these 
data gives an estimate of the Luttinger parameter Kc vs U/t (using 
a log scale starting at U/t = —2). 
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FIG. 17: (Color online) Bond kinetic energy modulation at the center 
of a chain of length L for various V at fixed U/t = 3. 



In the critical phase, the von Neumann block entropy gives 
access to the central charge c, and is consistent with a c = 1 
Luttinger liquid as expected (data not shown). In the gapped 
SP phase, the EE scaling is consistent with a saturation for 
large blocks. 

The transition from the critical BCS to the SP phase can be 
located when 'P{x) becomes exponential, or by looking at the 



bond kinetic modulation scaling. In Fig. 17 we plot the bond 



kinetic energy difference at the center of the chain as a func- 
tion of the chain length L. We can clearly see a finite value in 
the SP phase for V/t = — 1.5 and U/t = 3 for instance, while 
our data are compatible with an algebraic power-law with ex- 
ponent 0.66 for larger \V\. Locating precisely the transition 
is difficult since we expect a BKT behavior at the transition, 
and a weakly opening gap when entering the SP phase; this 
results in some uncertainty on this transition line in the phase 
diagram. Using this procedure, we have determined approx- 
imately the phase transition line shown in Fig. 13 Although 



our data are not very accurate, our numerical findings are in 
agreement with the low-energy approach: the BCS to SP tran- 
sition occurs for a. finite negative V for fixed [/ > 0. 

Note that the low-energy approach predicts that the tran- 
sition occurs when Kc reaches 1/3. According to our fitting 



procedure (see Fig. 16 1, this gives a similar estimate for its lo- 
cation. According to this value, the bond kinetic energy mod- 
ulation should scale as \/\fL at the transition, while we have 
measured a different exponent. In fact, it is known that loga- 
rithmic corrections are expected at this transition, and indeed 
our data can as well be fitted with a 1/^/L log L law. 

Moreover, along this transition line and from the low- 
energy approach, we expect an emergent SU(2) symmetry 

(SU(2)^) that should be reflected in identical exponents for 
V{x) and M{x). Fig. 14 b) displays our data in this region. 



and we do confirm a good agreement between the two expo- 
nents (compatible with Kc ~ 1/3). 

Concerning the transition from the critical BCS to the CDW 
phase, our data are compatible with a gap opening as soon as 
V > NU, in perfect agreement with the low-energy predic- 



tion (see Sec. IIID3i. Finally, for the same reasons as in 
the N = 2 case (see Sec. |IV C| l, we could not investigate the 
nature of the quantum phase transition between SP and CDW 
phases. We found that this transition is located in the vicin- 
ity of the V — — 9[//2 line (see Fig. 13 1. Unfortunately, as 
akeady stressed in Sec. II, we were not able to determine the 
symmetry contents of this line. 

As a final remark about the BCS phase, while quasi-edge 
states can be observed along the SU(2)c line or close to it (see 
previous section), they no longer exist deep in the BCS phase 
(for instance U/t ~ and V/t = — 6, data not shown). This 
might be understood from the strong-coupling regime using 
the mapping to a spin-3/2 chain with single-ion anisotropy: 
for large enough D > 0, the relevant low-energy states consist 
in Sf = ±1/2 on each site, thus leading to an effective spin- 
1/2 chain in its critical phase. In this region, we do not expect 
any edge physics as is observed numerically. We have not 
investigated in details the crossover between both regimes, but 
it could be easily answered by studying directly a spin-3/2 
anisotropic chain. 



VI. PHASE DIAGRAM IN THE iV = 4 CASE 

In this section, we investigate the phase diagram of 
model ([T) when A^ = 4 and in the {U/t, V/t) plane. 

From a technical point of view, since the local Hilbert space 
is quite large, we map the one-dimensional model onto a gen- 
eralized 4-leg Hubbard ladder with generalized rung interac- 
tions. Thus, we reduce the complexity of the DMRG algo- 
rithm, but we have to use a ID path going along the ladder. 
We have checked that the symmetry between the chains is al- 
ways restored during the simulations. Typically, we keep 2000 
states in our simulations and use OBC. 
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FIG. 18: Numerical phase diagram obtained by DMRG in the N = A 
case with L = 30. 



Fig. 18 shows the phase diagram for A^ = 4, obtained on 



a i = 30 chain. As expected, five phases are present: on the 
SU(2)c line in the attractive part of the phase diagram, there is 
the HI phase and, close to it, the critical BCS one arises. As 



expected from the strong-coupling argument (see Eq. ( 15 i), it 
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is followed by the RS phase and, on the other side of the HI 
phase, a CDW phase is stabilized. On the repulsive side, we 
detect a SP phase that was predicted in Sec. Ill C. We observe 
a good agreement with the low-energy prediction (see Fig.lTll. 



A. Properties along the SU(2)c line 

We start by looking at the SU(2)c hne V ^ 411 on the at- 
tractive side. As expected from the strong-coupling argument 
and from the low-energy analysis, the model should behave 
as an effective antiferromagnetic spin-2 Heisenberg chain, i.e. 
be in a Haldane phase. 

We have some evidence for such a HI phase thanks to the 
presence of (charge) edge states when OBC are used. Con- 
cerning the charge gap, in order to get the bulk result (and 
avoid edge states effect), one needs to compute Eo{Nf ~ 
4L + 6) - Eo{Nf = 4L + 4). Extrapolating our data on 
L = 16 and L = 32 chains for U = -2t and V = -8t, 
we obtain an estimate of 0.0038i, which is extremely small. 
Nevertheless, using the strong-couphng expression of the ef- 
fective exchange J^s = l/18t and the known Haldane gap*** 
A ~ 0.089 Jeff ~ 0.0049i, we get a finite gap of the same 
magnitude. 

Moreover, the pairing correlations (which correspond to the 
transverse spin correlations in the spin language) shown in 



Fig. 19 exhibit a short-range behavior compatible with a finite 



correlation length and a finite gap. 

However, since the correlation length of the spin-2 chain is 
known to be very large^^l (^ ~ 50), we will not try to char- 
acterize further this HI phase (by measuring its string order 
for instance), but the strong-coupling argument ensures that 
HI phase exists in some finite region of the phase diagram 
around the SU(2)c line. 



B. Critical BCS phase 



In Sec. 



II B we have argued why for fixed V/t, increas- 



ing U/t gives an effective single-ion anisotropy denoted D. 
For the spin-2 chain, it is known*^ that such a D term leads 
to an extended critical XY phase with central charge c = 1 
for 0.04 < D/J^ff < 2.4. Using our strong-coupling esti- 
mate and assuming a fixed effective Jeff, this would predict 
an extended XY phase for -1.996 < U/t < -1.82 for fixed 
V/t = -8. 



Fig. 19 shows numerical data for the pairing correlations 
obtained for a chain with L = 64. Indeed, algebraic behav- 
ior is observed close to the Haldane phase, in a region corre- 
sponding to -1.99 ^ U/t ^ -1.9 for V = -8t, in agree- 
ment with our strong-coupling estimate. 

In order to be more quantitative about this Luttinger liquid 
phase, and make a connection with the low-energy analysis. 








(a) - 


\^^^^ 




\X^^5*H^ 




\ ■« V^**^ 




\ V ' ^"^^I^ 


'SSih^ 


^\ > ',^^ 


'^C'SSi. - 


W^" 


^"x^ 


\ \ 


x^^ 




< — ,U/t = -2 


s X ^ 




—.[//( = -1.99 


\ "- 




^^l7/f=-1.98 






^-r//t=-1.96 


'% 




^^U/t = -l.'M 


\ ■ - 


■ 


u/i=-im 


\ 


- 


. .(//( = -L9 


\ - 







0.01- 




FIG. 19: Pairing correlations for A'^ = 4 at fixed V/t — —8 and 
L = 64. (a) Critical correlations are observed close to the HI phase 
represented by U/t — —2 on the SU(2)c line, (b) Both in HI and in 
RS phase, the pairing correlations are short-ranged. 



we use the asymptotics from Eq. (70 1 



V{x) 



{PUL/2 + x)PooiL/2))^x 
{n{L/2 + x)n{L/2)) - {n{L/2 

-^ + (-irAexp(-x/0, 



l/NKa 

x)){n{L/2)) 
(78) 



and use it to extract the behavior of Kc in the c = 1 gap- 
less phase. Note that we have measured the correlations 
from the center of the chain in order to minimize size effects 
due to OBC, and in the critical phase, we plot our data vs 

x' — d{x\L + 1)/ yJcos{-Kx/{L + 1)) in order to be able to fit 
over the whole range. '^^ 



Density correlations are shown in Fig. 20 a) in the critical 
phase, and they exhibit an algebraic decay. The anomalies 
are due to the subleading short-range staggered contributions 
since we plot the absolute values of A/'(a;), so that density cor- 
relations become difficult to fit close to the HI phase. We 
have fitted both correlations using the expressions above in 



Eqs. (78 1 and the resulting Luttinger parameter Kc is plotted 
in Fig. 20 b). We have an overall good agreement between 



the two independent fits, and we confirm the expected behav- 
ior that, starting from the HI phase, K^. first increases rapidly, 
and then diminishes when U/t increases. The gapless phase is 
characterized by Kc > 1, which corresponds for V/t = — 8 to 
-1.999 < U/t < -1.93. Thus, we conclude, from Eqs. (|78]l. 



that the leading instability is the BCS singlet-pairing. The ex- 
tended gapless phase is thus a BCS phase which differs from 
the one in the iV = 3 case by the fact that the staggered part 
of the density is short-ranged. 

For the choice of V/t = —8, we see that both pairing and 
density correlations become short-range when U/t ^ —1.9, 
where a RS (large-D) phase starts. As can be seen from the 
behavior of Kc, this corresponds to Kc = 1 = 4/A^ which 
is the criterion for the opening of the gap in the bosonization 



analysis done in Sec. Ill D 4 
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C. Quantum phase transitions 

Once we have determined the five phases in the phase dia- 
gram, we would like to clarify the nature of the quantum phase 
transitions. 

Starting from the SU(2)c line with V < and decreasing U 
leads to a CDW phase. According to the low-energy approach, 
and as confirmed numerically for A^ = 2, we expect a c = 
1/2 second-order Ising phase transition. However, due to the 
large correlation length in the HI phase, we cannot get reliable 
results. It would be much easier to check this criticality, as 
well as locate the critical D/J, by studying directly the spin- 
2 chain with single-ion anisotropy. 

On the other side of HI phase, i.e. increasing U, the phase 
transition to BCS critical phase was predicted to be in the BKT 
universality class with c — 1. Scaling of EE for L = 64 with 
V/t = -4.0 and U/t = -0.9 leads to c = 0.93 (when keep- 
ing m = 4000 states in the DMRG simulation). From the 
BCS phase and increasing U/t at fixed V/t < 0, our correla- 
tion functions in Figs. T9]|20]are compatible with a transition 
to a fully gapped RS phase when Kc becomes smaller than 1 
as expected. 

The transition from RS to SP is difficult to characterize due 
to finite-size oscillations in the quantities (including EE), but 
we have determined that it is located at a. finite negative V for 
fixed [/ > as found in the low-energy section (see Fig.lTll. A 
similar conclusion can be made for the transition from SP to 
CDW, which is located in the opposite quadrant as expected. 



VII. CONCLUSION 

We have established the zero temperature phase diagram 
of multicomponent (2A^-component) fermionic cold atoms, 
loaded in a ID optical lattice, at half filling. This entire work 
was done under the hypothesis that only contact interactions 
matter and that the interactions channels can be reduced to 
two: one singlet channel and one non-zero spin channel. The 
former hypothesis is very reasonable in the context of opti- 
cal lattices, and could be relaxed without affecting our main 



conclusions. The latter hypothesis requires N — 2 indepen- 
dent fine-tunings, and is therefore quite restrictive for large 
N. However, it requires no fine-tuning for N < 2, and should 
not be out of reach for moderate A^ < 4. 

As soon as iV > 1, we found that the phase diagram has a 
rich structure due to the degeneracy of the atomic states and 
the absence of spin-charge separation at half-filling. Several 
nonequivalent Mott-insulating phases emerge. Two phases are 
present irrespective of the value of A^ > 1: the SP and CDW 
phases, that both break translational invariance and are two- 
fold degenerate. We exhibited a hidden pseudo-spin SU(2) 
structure, involving spin-singlet, charged degrees of freedom, 
that generalizes a structure noticed long a go for the N — 1 
case in the context of the Hubbard model-I^Sl When special- 
ized to one space dimension, this structure yields a Haldane 
conjecture for attractive interactions: we show that such a 
system realizes a Heisenberg antiferromagnet of magnitude 
S = N/2, and, as a consequence, displays an alternating 
gapped (insulating) / critical (BCS superfluid) behavior ac- 
cording to the even / odd parity of N. We have found that this 
parity effect has an influence on large portions of the phase di- 
agram, and that ultracold fermions with N even can disclose 
two more insulating states: the HI and RS phases which are 
non-degenerate and display non-local string orderings. The 
A^ = 1 , 2 cases turn out not to be the generic cases of the 
odd/even families. Precisely, whereas for N — 2, only a crit- 
ical quantum phase transition occurs between the two non- 
degenerate insulating phases, for even A^ > 2, an intermediate 
gapless BCS phase arises between them. 

On top of the even-odd scenario, and within the low-energy 
approach, we found a subtle effect depending on the parity of 
A^/2. When A^/2 is odd, the HI and RS phases correspond to 
different phases and can be distinguished by string-order pa- 
rameters. In particular, the HI phase with odd N/2, i.e. odd 
spin, is an exemple of a topological ordered phase with N/A 
edge states. In contrast, when A^/2 is even, the HI and RS 
phases are related at low-energy by a duality symmetry and 
share the same ground-state properties. In this respect, the HI 
phase with A^/2 even, i.e. for even spin, is not topologically 
protected by its edge-state structure but is equivalent to a topo- 
logical trivial insulating phase, i.e. the RS phase. Thus, within 
the low-energy approach presented in this paper, our findings 
confirm the recent conjecture of Ref|33] 

In the light of the recent experimental achievements where 
cold fermionic gases with several components could be stabi- 
lized as highly symmetry systems, ^^ we hope that it will be 
possible in the future to unveil part of the richness that we 
highlighted in this work. In particular, the disclosure of the 
HI phase would be extremely important, as it displays exotic 
characteristics that have attracted a lot of attention in the past 
years and still does nowadays. 
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